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We study a periodic medium driven over a random or periodic substrate, characterizing the non- 
equilibrium phases which occur by dynamic order parameters and their correlations. Starting with a 
microscopic lattice Hamiltonian, we perform a careful coarse-graining procedure to derive continuum 
hydrodynamic equations of motion in the laboratory frame. This procedure induces nonequilibrium 
effects (e.g. convective terms, KPZ nonlinearities, and non-conservative forces) which cannot be 
derived by a naive Galileian boost. Rather than attempting a general analysis of these equations of 
motion, we argue that in the random case instabilities will always destroy the LRO of the lattice. 
We suggest that the only periodicity that can survive in the driven state is that of a transverse 
smectic, with ordering wavevector perpendicular to the direction of motion. This conjecture is sup- 
ported by an analysis of the linearized equations of motion showing that the induced nonequilibrium 
component of the force leads to displacements parallel to the mean velocity that diverge with the 
system size. In two dimensions, this divergence is extremely strong and can drive a melting of the 
crystal along the direction of motion. The resulting driven smectic phase should also occur in three 
dimensions at intermediate driving. It consists of a periodic array of flowing liquid channels, with 
transverse displacements and density ("permeation mode") as hydrodynamic variables. We study 
the hydrodynamics of the driven smectic within the dynamic functional renormalization group in 
two and three dimensions. The finite temperature behavior is much less glassy than in equilibrium, 
owing to a disorder-driven effective "heating" (allowed by the absence of the fluctuation-dissipation 
theorem). This, in conjunction with the permeation mode, leads to a fundamentally analytic trans- 
verse response for T > 0. 

PACS: 64.60. Fr, 74.20.D 



I. INTRODUCTION 



Non-equilibrium driven solids and liquids arise in a 
wide variety of different physical contexts. A common 
means of driving is to apply a constant or low frequency 
spatially uniform shear (either a constant shear rate or 
shear stress), which has beenjStudied extensively in col- 
loidal and polymeric systemsld. Such a driving is in a 
sense severe, since it is incompatible with a macroscopi- 
cally ordered solid, requiring the continual breaking of a 
non-zero density of bonds per unit time. If translational 
symmetry is broken explicitly by the presence of, e.g. a 
periodic substrate or quenched impurities fixed in space, 
a gentler sort of driving is possible. In this case, even 
a uniform translation of the system is nontrivial, and it 
can be driven out of equilibrium simply by applying a 
constant force or pulling at a constant velocity. 

A considerable number of such systems have been sub- 
jects of recent investigation. JHiese include flux lat- 
tices in type-II superconductorsoj, charge density waves 
(CDWs) in anisotropic conductorsQ, magnetic bubble 
array^l, and the magnetically-induced Wigner crystal in 



a two-dimensional electron gagj'Ll. In all these systems 
the relevant degrees of freedom - be they vortices in su- 
perconductors or electrons in metals - form a lattice in- 
side a solid matrix, provided by the superconducting or 
conducting material. Both a periodic potential (due to 
the underlying crystal lattice) and a quenched random 
one (due to material impurities and defects) are gener- 
ally present, though their relative importance can vary 
from system to system. Closely related problems also 
arise in microscopic models of friction and lubrication, in 
which a surface or monolayer is brought into contact with 
another surface and forced to slide |ijelative to it. Some 
recent simple models of earthquakeso, in which two elas- 
tic half-spaces are slowly driven past each other, also fall 
into the general class of driven disordered periodic elastic 
systems. 

Much of the recent focus in these pinned elastic me- 
dia has been on equilibrium behavior, since, unlike their 
colloidal and polymeric counterparts studied in the shear 
geometry, these systems exhibit a nontrivial competition 
between the external (substrate or disorder-induced) po- 
tential and the tendency for local order. In the random 
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case, this was argued by LarkinB to generate long-range 
elastic distortions. More recent works have reinvestigated 
this problem in some detail, suggesting the existence of 
a novel "Bijagg-|glass" phase in three and possibly two 
dimensions .tHTla In the periodic case (known as intrinsic 
pinning in the vortex community), the potential can ci- 
ther lock-in commensurate phases, generate finite (but 
qualitatively unimportant) incommensurate distortions 
of the lattice, or stabilize anisotropic liquid- crystalline 
states. Though for both types of matrix some detailed 
questions remain unanswered (in particular, the stabil- 
ity of the studied elastic glassy phases to proliferation of 
topological defects), the equilibrium phases and transi- 
tions in these systems have been extensively studied and 
are fairly well understood. 

Once the elastic medium is driven, however, a host 
of new questions arises: Under a uniform applied exter- 
nal force /ext., what is the mean velocity v{fcxt.) (the 
IV curve, in the context of superconductors and CDWs), 
and is this uniquely determined or dependent upon his- 
tory or other variables? Do ordered solid phases exist at 
low temperatures or weak disorder, and if so how are 
they characterized? Can one develop non-equilibrium 
phase diagrams, with phases classified by order param- 
eters and symmetry analogous to equilibrium problems? 
What are the properties of the resulting non-equilibrium 
dynamic phase transitions and the nature of the fluctua- 
tions? Under what conditions does such a system reach 
a non-equilibrium steady state, and what are the hydro- 
dynamic modes in this case? 

Along with these new questions come a number of new 
physical variables which play important roles. While 
equilibrium behavior is relatively insensitive to detailed 
dynamical laws, there is no reason to expect this for 
a driven system. In particular, the nature of dissipa- 
tion surely plays a crucial role. This distinguishes, e.g. 
friction and lubrication, in which (to a good approxima- 
tion) energy can only be redistributed among vibrational 
modes of the solid, from charge-density-wave or vortex 
solids, in which energy can be transfered out of the collec- 
tive modes to dissipate into electronic degrees of freedom. 
In the latter case, which we focus on here, it is appropri- 
ate to consider overdamped dynamics, while for the for- 
mer inertial effects may be significant. Also crucial are 
conservation laws, which generally give rise to additional 
hydrodynamic modes. In equilibrium, the interactions of 
these with elastic degrees of freedom are constrained by 
the fluctuation-dissipation theorem (FDT) not to modify 
static correlations. For the nonequilibrium lattice, how- 
ever, they must be treated explicitly. Finally, a driven 
system may exhibit complex dynamics even in the ab- 
sence of external "thermal", or time-dependent, noise. 
When the external noise level is low, this deterministic 
dynamics can give rise to steady, periodic, quasiperiodic, 
or (spatio-temporally) chaotic solutions. While in the lat- 
ter case, dynamical chaos likely gives rise to an effective 
"temperature" and restoration of ergodicity, more regu- 
lar solutions need not explore the full phase space and 



can in principle behave with almost arbitrary complex- 
ity! Of course, factors such as the nature of the external 
potential, dimensionality, and range of inter-particle in- 
teractions, which control the equilibrium state, influence 
the driven dynamics as well. 

A complete answer to these questions for all possi- 
ble cases is beyond the scope of a single paper (or re- 
searcher!). Instead, we will focus here on systems with 
overdamped, dissipative dynamics which reach statisti- 
cally steady states. This can occur either due to the pres- 
ence of external time-dependent noise, which forces the 
system to explore the available phase space, or due to in- 
trinsic chaotic dynamics. Note that this condition is vio- 
lated by pY^rdamped phase-only models (e.g. Fukuyama- 
Lee-Rical3) of CDWs, in which there is a uniaue periodic 
long-time attractor, as shown by MiddletonJiil It may also 
be violated in zero-temperature simulations of vortex dy- 
namics under some circumstances, as observed recently 
by Nori ip. some regimes for a strong periodic pinning 
potential. EJ We expect, however, that in the presence of 
a random potential chaotic dynamics is much more ger- 
mane, and furthermore, that most systems of experimen- 
tal interest contain appreciable external thermal noise. 
Nevertheless, where possible, we will comment on the 
extensions of our conclusions to the noiseless case. 

Our approach to the problem is first to devise a means 
of classifying the non-equilibrium phases, then to deter- 
mine the dynamical equations of motion governing these 
phases. Finally, using these equations, we can begin to 
calculate the phase diagram and the properties thereof. 
We summarize the main results below, reserving a com- 
parison with previous work for the discussion section of 
the paper. Some preliminary versions of some of the re- 
sults of this paper appeared in Ref. ^ 



A. Order parameters and correlation functions 

A framework for classifying phases was introduced in 
Ref. |l^ in the context of CDWs, and we generalize it 
here to more complex periodic media. As in equilibrium, 
phases are characterized by broken symmetry. An ordi- 
nary crystal is globally periodic, and hence has broken 
translational and rotational symmetry. Because the pe- 
riodicity is sustained over long distances, we say that 
it has translational long range order (LRO). To quantify 
this notion, we must introduce suitable order parameters. 

We consider initially a general model where the con- 
stituents of the driven lattice may be "oriented mani- 
folds" , which are extended in some direction, such as 
magnetic flux lines in a three-dimensional superconduc- 
tor. The configuration of such a lattice is described by 
labeling each manifold by its undisplaced equilibrium 
transverse position, x. The internal coordinates parallel 
to the oriented manifold (e.g., along the magnetic field 
direction in the vortex lattice) are parameterized by a 
(i;-dimensional vector z, and r = (x, z). The number of 
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transverse dimensions is denoted by dt and the dimen- 
sionality of space \s d = dt + di. The lattice may be 
described by a density field smoothed out on the scale of 
the lattice spacing, /9(x, z) , which in the ordered phase 
becomes periodic in the transverse direction. For the vor- 
tex lattice the density is the local magnetic induction, 
while for the CDW and the Wigner crystal it is the elec- 
tronic charge density. In the ordered phase the collection 
of oriented manifolds acquires long-range periodicity, de- 
fined by a reciprocal lattice with basis vectors {Q}. The 
density field is then written as 

p(x,z)=po + ^5]pQ(r)e^^■^ (1.1) 
Q 

where po is the mean density. The complex Fourier 
components PQ(r) satisfy P_q = pq. All the ampli- 
tudes transform aS|-AQ — > pqe**^ '* under a lattice trans- 
lation X ^ X -I- atZl. A nonzero expectation value of 
the pq's thereby indicates broken translational symme- 
try. As a result of the broken translational symmetry, 
long-wavelength fluctuations in the ordered phase can be 
described in terms of long-wavelength distortions of the 
Fourier modes, which can in general be written in terms 
of an amplitude and a phase as 

PQ(r)=PQ,o(r)e^«-*W. (1.2) 

Note that the phase in each Fourier amplitude is not an 
independently fluctuating variable in the ordered state. 
It is constrained by, e.g. cubic interactions of the pq, 
which leave only the vector of phases </>. They corre- 
spond to distortions of the lattice, and indeed <j) can be 
interpreted as a sort of displacement field (but see below). 

So far, this discussion applies equally well to both 
equilibrium and driven solids. Consider now the time- 
dependence of the local density. In a moving solid, we 
expect cf) = vt + (f), so that 

p(x, z, i) = po + ^ E '^q(^' i)e*'*-[^+^*l, (1.3) 
Q 

where 

PQ(r,t)=pq(r,i)e-^Q-*. (1.4) 

Physically, the oscillations in the density simply reflect 
the fact that individual constituents of the lattice pass 
any given point in a regular periodic fashion. The set 
of Pq fields comprises the order parameters for the non- 
equilibrium system. Neglecting topological defects (i.e. 
dislocations) amounts to assuming that the amplitude 
Pq,o is constant. Elastic deformations of the solid are 
then described entirely by phase fluctuations in terms of a 
single- valued coarse-grained displacement field 0(r). It is 
important to distinguish between the displacement field 
defined in the moving or crystal frame (where the con- 
ventional phonon field is defined) and the displacement 



field in the laboratory frame. Throughout this paper we 
denote by </>(r) the displacement field in the laboratory 
frame, while we reserve the symbol u(r) for the conven- 
tional phonon field. Many of the experimental systems 
of interest form triangular lattices. In this case, (pq) ^ 
for three shortest wavevectors in the ordered phase (the 
brackets indicate a time average or average over thermal 
noise) . 

A quantity of immediate experimental interest is the 
static structure factor, which is the Fourier transform of 
the equal time density-density correlation function and 
is given by 

5(q) = (IMq,i)P> 

= J2 [ e^(''-'^)--(pQ(r,t)pJj(0,i)). (1.5) 

Q 

Here q — {q_t,<iz) is the full wavevector. The order pa- 
rameter correlator, 

CQ(r,t)- [(pq(r,i)p^(O,0))]_ , (1.6) 

governs the behavior of the structure factor. We used 
angular (square) brackets to denote thermal (disorder) 
averages. If CQ(r,0) const, for large \r\, S{q) has 
a sharp (resolution limited) Bragg peak a,t qt — Q- If 
CQ(r,0) ^ these broaden to power-law form. If 

CQ(r,0) decays exponentially for large |r|, only diffuse 
peaks if any exist. In these three cases, the system is 
said to possess long-range, quasi-long-range, and short- 
range translational order (at wavevector Q), respectively. 
These will be abbreviated LRO, QLRO, and SRO, respec- 
tively, in the remainder of the paper. 

Cq (r, t) also describes the extent of temporal peri- 
odicity. By analogy to the translational characteriza- 
tion, we say that the system possesses long-range, quasi- 
long-range, or short-range temporal order (at frequency 
w = Q . v) if Cq(0, i) goes to a constant value, decays 
algebraically, or decays exponentially, respectively. This 
temporal order can be experimentally probed through 
the dynamic structure function 

S{q,L0)^{\Sp{q,u)\') (1.7) 
= ^ /■ e^(i-Q)-'(-Qv)t^^Q(r,^)^^(o,0)). 

Q -^r,* 

Systems with temporal LRO will display sharp peaks at 
multiples of the washboard frequency w = Q . v in the 
5(q, w) and in the related power spectrum of velocity 
fluctuations. Such system should also exhibit "complete 
mode locking" to an arbitrarily weak external periodic 
driving at these frequencies. 

In a random system, temporal order is generally more 
robust than the translational order. Physically, the dif- 
ference arises because impurities inhomogeneously stress 
the system. The responding non-uniform distortion, 
however, can have very little fluctuation in time, and 
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thereby can leave the temporal ordered relatively unaf- 
fected. It will, of course, have some effect on the dynam- 
ics, because disturbances propagate differently atop the 
non-uniform background, and because the local strains 
lower the barriers to defect nucleation. We cannot ex- 
clude the possibility of a phase in which dislocations 
are unbound, but frozen in the laboratory frame. Such 
a phase would exhibit translational SRO but temporal 
LRO, and would be the driven analog of the vortex glps 
phase originally proposed by Fisher, Fisher and Huseli3. 
Given the relative instability of this phase in equilibrium, 
we think this scenario is, however, somewhat unlikely. 
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Fig.l:(a) Schematic illustration of a moving glass, which 
becomes unstable to the proliferation of dislocations 
with Burgers vectors along the direction of motion. At 
large velocities, we expect that unbound dislocations 
will be widely separated, with spacings Xsm and ysm m 
the directions longitudinal and transverse to the motion, 
respectively. A naive estimate based on the linearized 
theory probably gives an upper bound for these scales: 
Xs-m ^ ceea^i-fvy/A'^ and ^ ceea^ {■-/v f / A'^ . Non- 
linear (KPZ) effects very likely shorten these lengths 
considerably. On larger scales, the "lattice" crosses over 
to the driven smectic state displayed in (b), with only 
transverse periodicity and liquid-like order along the di- 
rection of motion. The dislocations with Burgers vectors 
along the direction of motion allow the weakly coupled 
liquid channels to move at different velocities. 

What particular types of translational order can in 
principle arise in a driven system? Of course, a disor- 
dered liquid state is possible (and may be the only stable 
phase in low dimensions). When translational order is 
present, it can occur at a variety of wavevectors Q. For 



weak disorder at low temperatures, it is natural to expect 
that a full reciprocal lattice of wavevectors characteriz- 
ing a crystal should be important (i.e. have appreciable 
Cq). In two dimensions, for instance, the smallest of 
these would typically be arranged into a hexagon. If, as 
seems natural, and can be shown in some simple models, 
reflection symmetry perpendicular to the applied force is 
not broken, this can take one of two orientations, with 
a diagonal oriented either parallel or perpendicular to 
the force. Since the applied force breaks rotational sym- 
metry, however, there is no reason for the correlations 
at all six points to be identical. Instead, if the system 
evolves continuously with increasing temperature or dis- 
order, it is natural to expect that LRO will be lost first 
at some subset of these wavevectors. The surviving state 
has a lesser periodicity, with only a single line of Bragg 
peaks (symmetry requires that the other solid peaks dis- 
appear pairwise) . It represents therefore a layered liquid 
or smectic state. This can be either a longitudinal smec- 
tic, with ordering wavevector parallel to the velocity, or a 
transverse smectic, with periodicity in the direction per- 
pendicular to the velocity. We argue below that only the 
case of a transverse smectic shown in Figure 1(b) is sta- 
ble, and indeed, that this is likely to be a more generic 
state than the true solid when the system is driven. 



B. Hydrodynamics and elasticity 

There are at least two analytical approaches to calcu- 
lating CQ(r,t). One is to construct a density- functional 
or Landau-like theory for the order parameters pq. For 
the non-equilibrium driven system considered here, this 
would be a set of stochastic partial differential equations. 
This has the advantage of allowing large amplitude fluc- 
tuations, and hence including dislocations in a natural 
way. The disadvantage of this approach is its intractabil- 
ity. It is often difficult to recover relatively simple prop- 
erties in the ordered phase. A second approach is to 
assume a large degree of local order, so that the ampli- 
tude of the order parameters fluctuates very little. In 
this phase-only or elastic approximation is nat- 

ural, and only the cj) fields remain in the description. 
This method has the disadvantage of excluding topologi- 
cal defects, which must be reintroduced by hand to com- 
plete the description, as has been successfully done for 
equilibrium systems, such as for example Kosterlitz apd 
Thouless treatment of 2d superfluids and 2d meltingcJ. 
Determining the relevance or irrelevance of such defects 
upon the elastic description in disordered systems is a dif- 
ficult and unsolved problem. If, however, the phase-only 
approximation predicts only small displacements in the 
ordered phase, it can be expected to be a self-consistent 
approximation out to a rather long length scale. If, by 
contrast, large distortions are found, the original assump- 
tion of order is inconsistent, and we have discovered an 
instability of the solid phase. 
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To determine the nature and stability of the possible 
ordered non-equilibrium phases, this paper focuses on 
the elastic approach. Assuming constant amplitude of 
the order parameters, Cq can be written in terms of the 
displacement field correlation as 



One noteworthy feature of Eq. 1.11 is the proliferation 



Cq(i":^) = IPQ.oI 



„iQ-[</>(r,t)-</'(0,0)]\ 



(1.8) 



where we have pulled out the amplitude factor. If is a 
Gaussian random variable. 



CQ(r,t) = \pqM'^ exp{-QiQjBij{r,t)}, 



(1.9) 



where i and j denote Cartesian components and Bijiv) 
are the components of a mean-square displacement ten- 
sor, 

B,,{v,t) = - 0.(O,O)][0,(r,t) - 0, (0, 0)])]^„^. . 

(1.10) 

Generally, the displacements fluctuate in a non-Gaussian 
manner, so that Eq. 1.10 is not strictly correct. However, 



we expect its qualitative implications to hold. If Bij (r, t) 
goes to a constant at long distances (times) , then so does 
CQ(r,t). If Bij{r,t) grows algebraically, then CQ(r,t) 
has stretched exponential form (though not necessarily 
with the naive stretching exponent), suggesting that in 
a physical system topological defects might proliferate 
and lead to short-range order, and if Bij (r, t) grows log- 
arithmically, then Cq (r, t) decays as or slower than a 
power-law. 

To calculate Bij{r,t), we employ the analog of hydro- 
dynamic equations of motion. In general, the hydrody- 
namics of systems far from equilibrium is far more com- 
plex than its equilibrium counterpart. In particular, fluc- 
tuations about the nonequilibrium steady state do not 
satisfy a fluctuation-dissipation theorem and the exter- 
nal driving force breaks both the rotational and reflection 
(parallel to the force) invariance of the equilibrium sys- 
tem. As a result, the hydrodynamic equations contain 
many more parameters than in equilibrium. A general 
construction based only on symmetry constraints is thus 
not very useful, and a concrete derivation, which provides 
precise values for these parameters, is desirable. We per- 
form such a derivation in Section ||. Our first main result 
is the complete non-equilibrium hydrodynamic equation 
of motion for the driven lattice, 

7y9t0j = Aiajda(l>t + Bia/3jdadf3(l>j (l-H) 

where the coefficients A,B, and C are non-zero only when 
the number of subscripts taking values of axes perpendic- 
ular to the driving force is even. Explicit expressions 
for these are given in Section ||. F and r] represent 
quenched random and external "thermal" forces, respec- 
tively. Eq. 1.11 remedies the deficiencies of various-ad- 



of gradient terms beyond the usual equilibrium elastic 
ones (contained in the B term). These represent convec- 
tive effects and dependence of the substrate-lattice inter- 
action upon the local deformation of the lattice. To prop- 
erly account for them, it is crucial to treat the phonon 
modes near the zone boundary, not considered in previ- 
ous calculations. Especially important are the A terms, 
which lead to propagating modes at low frequencies and 
wavevectors. 

Another key feature is the random force Fi[r, 0, t]. It 
is distinguished from the equilibrium form in two ways. 
First, it contains non-trivial time-dependence, as can be 
seen from 



Fdr,(t>,t] 



(1.12) 



where Fi (r) is a time- independent quenched random vari- 
able, which we will refer t o as a static pinning force (the 
exponential factor in Eq. 1.12 comes from the transfor- 
mation from the moving to the laboratory frame). The 
terms with Q-v ^ therefore oscillate in the sliding solid. 
The second distinction is seen from the decomposition 



F,ir)^Frir)+Frir). 



(1.13) 



The first term on the right hand side of Eq. 1.13| is the 
equilibrium component of the static pinning force and can 
be written as the gradient of a potential, as required by 
the fluctuation-dissipation theorem. Its correlations are 
approximately given by 



[FP{r)F';^{r% 



-5,5jr(r-r'), 



(1.14) 



where r(r) is the correlator of the Gaussian random po- 
tential. The second term in Eq. 1.13 is the nonequilibrium 
part of the static pinning force, with correlator 



[Frir)Fpir')]^ 



(1.15) 



The variance gij of this static force is given in Eq. 2.51 
and vanishes in the absence of external drive. A non- 
zero gij implies that F^°'^ is non-conservative, violating 
the fluctuation-dissipation theorem. 



C. Analysis 



hoc equations of motion proposed in previous works 



A general analysis of Eq. 1.11 is quite difficult. In 
principle, the stability and behavior of the putative 
moving glass can be determined by an RG analysis of 
the full equations of motion (Eqs. 1.11), including both 
the KPZ nonlinearities (Ciajpk) and the random forces 
{Fi{r, 4>,t)). Past ejsqjiexience with other nonequilibrium 
KPZ-like equationsll3G3 suggests that instabilities are 
quite ubiquitous in low dimensions. We expect such in- 
stabilities will occur also in this case, at least in two di- 
mensions and quite possibly in three, and leave an analyt- 
ical check of this hypothesis to younger researchers who 
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may have enough career years ahead of them to complete 
the calculation. Provided such an instability occurs, can 
any residual order survive in the driven state? Two phys- 
ical realizations of partially ordered moving states have 
already been suggested: the longitudinal and the trans- 
verse smectic. A longitudinal smectic is equivalent to 
a conventional driven CDW, studied earlier by Chen et 
al.E3 and by Balents and FisherE^I. These authors con- 
cluded that this phase is unstable in two dimensions, 
but may exist at large velocities in three dimensions (al- 
though the role of dislocations and KPZ nonlinearitics in 
three dimensions deserves further study). 
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Fig. 2: (a) The static structure factor for the moving glass 
sketched in Fig. 1(a), with sharp peaks at all the recip- 
rocal lattice vectors of a triangular lattice. The static 
structure factor of the driven transverse smectic is shown 
in (b). Sharp peaks (dark spots) only appear along the 
Qy axis; some remnant of short-range local triangular 
lattice order can appear in the form of weak and diffuse 
peaks at finite (light spots). 

The only possible ordered phase in d = 2 is thus a 
transverse smectic (discussed in the previous subsection 
and illustrated in Fig. 1(b)) with some degree of order at 
a periodic set of Q perpendicular to the velocity. This is 
of course also a possibility in d = 3, regardless of the sta- 
bility of the lattice. The marginal stability of the driven 
lattice in three dimensions allows for a dynamic phase 
transition between the driven smectic (at intermediate 
velocities) to a moving lattice (at high velocities). We 
illustrate the corresponding dynamic phase diagrams in 
Fig. 3. 

The latter part of the paper is devoted to an analy- 



sis of this possibility. In Section we present the hy- 
drodynamic equations of motion for the smectic. These 
include a simplified version of Eq. |l . 1 1| , supplemented 
by an additional one for the conserved particle density. 
Unlike in a vacancy/interstitial-free solid, this is not gen- 
erally slaved to the smectic displacement, and constitutes 
a separate hydrodynamic modeS This is the well-known 
permeation mode in smectic liquid crystals. 

Even these equations are somewhat intractable, so in 
Sections |v|-[v| we study the "toy model" in which the 
permeation mode is decoupled from the smectic displace- 
ment (jjy = (f'- This is best done using renormalization 
group (RG) techniques, which are controlled in two lim- 
its. At and near d = 3, the RG is controlled by a low- 
temperature fixed point, which is analyzed using a func- 
tional RG in Section ^ The fixed-point temperature T* 
increases with decreasing dimension until in d = 2, more 
conventional sine-Gordon RG techniques can be applied. 
Directly in three dimensions, there is a slow asymptotic 
approach to zero temperature. 

From these RG calculations come several concrete pre- 
dictions. In a three-dimensional smectic, the structure 
factor has power-law Bragg peaks along the axis in mo- 
mentum space perpendicular to the velocity, as illus- 
trated in Fig. 2(b). Because these peaks are entirely 
transverse, the smectic exhibits neither narrow-band- 
noise nor complete mode locking (although incomplete 
mode locking is of course possible for sizeable ac drives). 
The response to a force transverse to the mean direc- 
tion of motion is a superposition of two effects. First, the 
permeation mode provides a non-vanishing linear com- 
ponent of the velocity v'J"™' — /ipcrm.f_L- Secondly, the 
smectic responds in a non-linear manner, resembling a 
threshold at low temperatures but crossing over to linear 
response at very low forces for T ^ 0. 

In two dimensions, the RG study of the elastic model 
also predicts only short-range asymptotic transverse 
translational correlations and quasi-long-range temporal 
correlations. These calculations actually indicate that 
the system is outside the regime of validity of the elas- 
tic approximation. We expect that in d = 2 an eventual 
instability of even the smectic state to unbinding of trans- 
verse dislocations may occur at larger length scales. Nev- 
ertheless, the RG results should hold up to this length, 
and, as in three dimensions, also predict superimposed 
linear and non-linear transverse velocity responses. 

We conclude the Introduction with a summary of the 
remainder of the text. Section II describes the deriva- 
tion of hydrodynamic equations of motion for the driven 
lattice, which we study perturbatively in Section III. Sec- 
tions IV presents the equations of motion for the trans- 
verse smectic, which are analyzed using renormalization 
group techniques in Sections V and VI. The specific pre- 
dictions of the RG for correlation and response functions 
are given in section VII. Section VIII summarizes our 
results, makes comparisons with other work, and gives 
a synopsis of remaining open questions and future ap- 
plications of these ideas. Finally, six Appendices give 
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technical details unsuitable for the body of the paper. 
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Fig.3:Schematic phase diagram for a finite temperature 
driven solid in (a) d = 2 and (b) d = 3 dimensions. In (a) 
we have used a dashed line as a boundary between the 
moving liquid and the transverse smectic to emphasize 
that the smectic might be unstable in 2d and therefore 
asymptotically indistinguishable from a moving liquid. 
In (b) we have similarly used a dashed boundary be- 
tween the moving smectic and lattice, to emphasize that 
it is likely that the moving lattice is unstable even in 3d. 
In (b) the hatched "?" region indicates the interesting 
possibility of a moving Bragg glass in 3d, at low drives. 



II. DERIVATION OF HYDRODYNAMICS 

As discussed above, the goal of this paper is to study 
the non-equilibrium steady states that arise in driven pe- 
riodic media. We will focus on the low-energy and long- 
wavelength properties of these steady states, in cases 
where these are uniquely defined. This should be the 
case provided ergodicity is achieved, either through a 



small non-zero "thermal" noise or via intrinsic chaotic 
dynamics of the system in the absence of external noise. 
We will not discuss zero-temperature systems in which 
relatively simple global attractors exist with non-chaotic 
(e.g. periodic) dynamics. 

In the limit of interest, then, we expect a sort of hy- 
drodynamic description to hold. Such hydrodynamics is 
particularly successful in equilibrium because it is highly 
constrained. It must respect both detailed balance (to re- 
duce to equilibrium statistical mechanics) and the sym- 
metries of the system - in this case translations, rota- 
tions, and refiections. Out of equilibrium, a putative 
steady state equation of motion can be much more gen- 
eral. The fluctuations around the driven state need not 
satisfy any fluctuation-dissipation relation, and the ex- 
ternal driving force breaks both rotational and reflection 
invariance. 

For this reason, there are many more parameters in 
the hydrodynamic description. In the absence of fur- 
ther input, it is therefore considerably less powerful than 
equilibrium hydrodynamics. To make it useful, we need a 
means of calculating these parameters "microscopically" . 
This is possible for weak external potentials. 

The calculation proceeds in two stages. Beginning 
with a microscopic lattice Hamiltonian, we first explicitly 
coarse-grain the eqiiations of motion. This is a mode- 
elimination procedure in momentum space, reminiscent 
of a single step of rescaling in a renormalization group cal- 
culation. It is this stage which requires a weak external 
potential, since the elimination can only be performed 
perturbatively in these mode-coupling terms. The end 
result is an effective equation of motion for the phonon 
modes with momenta within a small sphere of radius 
A <C 27r/a in momentum space (a is the lattice spacing). 

The second step is to transform the equation of mo- 
tion from the moving frame (in which the conventional 
phonon coordinate or displacement field is defined) to 
the laboratory frame (more relevant for physical mea- 
surements). While technically much simpler than the 
previous mode-elimination step, it is possible only at this 
stage, since the transformation can be carried out only as 
a gradient expansion, with the small parameter Aa<^l. 



A. Formulation 



The formulation of the problem begins with some mi- 
croscopic Hamiltonian, describing particles connected by 
springs. The natural degrees of freedom are thus the 
displacements of these particles, whose gradients are 
expected to be small provided the potential acting on 
the particles is weak. The high-momentum modes be- 
ing eliminated here therefore describe small relative dis- 
placements of nearby (say, neighboring) particles. Small 
though these are, they are crucial to the physics of the 
sliding solid. 
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Furthermore, the effects of the fixed external potential 
are expected to be most pronounced when it is has strong 
periodic components commensurate with the driven lat- 
tice. This is intuitively reasonable, and indeed comes out 
of our more detailed calculations. This physics, however, 
comes precisely from these high-momentum modes at the 
scale of the lattice spacing. Even if the lattice is not 
pinned into a static configuration, these modes are the 
most strongly coupled to the static matrix, and thereby 
give rise to the seeds of interesting non-linear dynamics. 

We consider a general model in which the constituents 
of the driven lattice may be "oriented manifolds" , which 
are extended in some directions. This allows for sys- 
tems including, for example, vortex lines in a three- 
dimensional superconductor. The conformation of a 
manifold with undisplaced equilibrium transverse posi- 
tion X in a distorted lattice is described by the displace- 
ment vector u(x, z), such that the true position is 



In Fourier space, using reflection symmetry, the elastic 
matrix is 



X(x, z) = x + u(x, z). 



(2.1) 



where u and x are dt-component vectors, while z has di 
components. 

A natural microscopic Hamiltonian, valid in most 
cases, is 

H = ^ y [x - x' -f u(x, z) - u(x', z'), z - z'] 

XX' 

+ ^ /'c/[x + u(x,z),z]. (2.2) 

X "'^ 

Here V[x, z] is a two-body manifold-manifold interaction, 
and C/[x, z] is the external potential of the static medium. 

Expanding the interaction potential in the usual way 
gives the elastic energy 

H^\Y. y" A-u(x-x',z-z')u'(x,z)7.J"(x',z') 

x,x' 

-f^ /c/[x + u(x,z),z], (2.3) 

X "'^ 

where the elastic matrix 

Kv(r) = Y. I d,d,V{v')6^,o5{7.) - d,d,V[v]. (2.4) 

X' 

At this point, it pays to establish some notation. In 
the previous equation, we have already adopted a conven- 
tion for indices. Latin indices alphabetically following i 
denote transverse (x) coordinates, while a,b, . . . ,h de- 
note longitudinal (z) ones. If an index may range over 
the full space, a Greek index a,(3, . . . will be used. We 
will attempt to use x and z (with primes, etc.) exclu- 
sively as transverse and longitudinal coordinate vectors, 
with corresponding momenta and . These will often 
be assembled into vectors in the full c?-dimensional space 
denoted r = (x, z) and q = (q*, qz). 



^ J Z ^ J Z 



(2.5) 



Overdamped relaxational dynamics is defined by 



^dtUi{Y,t) = - 



SH 



(2.6) 



where f is an external force. This is equivalent to the 
equation of motion 

X' 

+F,[x + u(r,<)+vi,z], (2.7) 

where we have shifted u to remove the external force, 
V = f /7, and 



F,(r) = -9.C/(r). 



(2.8) 



In general, v — f/7 is not the true steady-state veloc- 
ity for a given force f , since interactions with the impu- 
rities intermittently pin the lattice and thereby reduce 
the sliding speed. This finite renormalization (which is 
quantitatively small at large velocities) can, however, be 
neglected for the current purpose of studying the prop- 
erties of the steady state. This is analogous to ignoring 
the mass renormalization in field theory (or Tc shift in 
critical phenomena). Of course from the experimental 
point of view, the velocity curve v{f) (like the true Tc) 
is an important and measurable quantity. 



B. Mode elimination 

Having set up the dynamics, we are now in a position 
to perform the coarse-grainingr. We will use the Martin- 
Siggia-Rose (MSR) formalismEJ, which allows the mode- 
elimination to be performed by functional integration. 
The object of interest in the MSR method is the statis- 
tical weight 

VF[u,u]=e-^, (2.9) 
where the action S ^ Sq + Si, with 

X "'2;* L X' 

(2.10) 

and 

Si = -J2[ ii»(r,t)i^,[x + u(r,t) +vt,z]. (2.11) 

^ Jz.t 
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With this choice of weight, correlation functions are 
given by the functional integral 



{u,{r,t)uj{r',t')---) = J [du][du]u^{r,t)uj{r',t') 



(2.12) 



In this expression, the left-hand side should be inter- 
preted as simply the product of t he specified Ui evalu- 
ated along the solution of Eq. 2.7. The right-hand side 
is the result of functional integration over all functions 
u, u. Because of this identity, we will freely employ an- 
gular brackets in either context. We note that this equal- 
ity relies crucially on proper regularization of the equal- 
time correlators in the field theory. In particular, we will 
choose the causal convention 



{u,{r,t)uj{r',t))=0. 



(2.13) 



Correlation functions involving the field u have the physi- 
cal interpretation of response functions, as can be seen by 
simply differentiating a correlation function with respect 
to an applied force. This convention then simply implies 
that there is no instantaneous (and hence discontinuous) 
response to a perturbation. 

Within the MSR formalism, we can now readily carry 
out the mode elimination. The fields are first separated 
into two parts via 



u{r,t) = u<(r,f) + u>(r,i), 
u(r, t) = u<(r, t) + u>(r, t), 

Here the slow and fast fields are defined by 



^ici-r—iLjt 



u<(i",i) = / / u(q, 

|qt I <A J <\z^ 

u<(r,i) = / / u(q, 

|qt I <A J c{z^ 

u>(i".^) = / / u(q> 

J |qt I > A J <\z^ 

u>(r,i)^/ / u(q,(^)e-*'i-+^"*, 

|qt I >A J <\z^ 



— -iq-r+icji 



^iq-r—iujt 



(2.14) 

(2.15) 
(2.16) 
(2.17) 
(2.18) 



where we have adopted opposite Fourier sign conventions 
for the displacement and response fields, and all momen- 
tum integrations are restricted to within the Brillouin 
zone. Correlation functions of the slow fields describe 
all the long-wavelength (i.e. hydrodynamic and elastic) 
behaviors of the system, and can be obtained from the 
effective weight 



off.[U<,U< 



= J [du. 



(2.19) 



For simplicity, let us consider the model at zero tem- 
perature - i.e., in the absence of any external time- 
dependent noise. Because the potential couples the slow 
and fast modes, we will nevertheless obtain non-trivial 
renormalizations of the slow dynamics from the mode 
elimination. 



We proceed by inserting the decomposition in Eq. 2.14 
into the action S. Furthermore assuming that the po- 
tential U is weak, we may expand the exponential and 
perform the functional integrations over the fast modes 
order by order, re-exponentiating the resulting expres- 
sions, which then depend only upon the slow fields. The 
elastic part of the Hamiltonian is diagonal in momentum 
space, 



5*0 = / Uiiq,uj) 



ijLU + Kij{q,LL!) 



u,{q,iu), (2.20) 



so that upon decomposition the slow and fast fields are 
decoupled in this term: 



S'o[u, u] = S'o[u<,u<] + 5o[u>,u>] 
The effective potential is therefore 

5eff = 5'o[u<,u<]-ln(e~^A , 
\ / o> 



(2.21) 



(2.22) 



where the angular bracket with the subscripts indicates 
integration over the fast modes (>) with the additional 
weight factor e^'^" (0). 

From this point onward, the treatment differs for the 
periodic and random potential, so we divide the remain- 
der of this subsection into two parts. Each involves the 
conceptually straightforward perturbative calculation of 
the average in Eq. 2.22. This is somewhat tedious tech- 



nically, so details will be given in Appendices A and B. 



1. Periodic potential 

We can specify a periodic potential by the Fourier de- 
composition. 



C/(x)=5]e^Q-[/Q. 
Q 



(2.23) 



Note that we have taken the potential to be independent 
of the continuous longitudinal coordinates z. 



Inserting the above Fourier decomposition and evaluating the expectation value in Eq. 2.22| to second order in U 
gives corrections to 5*0 which correspond to a renormalized equation of motion (see Appendix A). Since for small A, 
the remaining "slow" fields have small gradients, the corrected equation of motion may be written in the continuum 
approximation. It takes the general form 



^ijdtUj = AiajdaUj + Biaf3jdadf3Uj + Ciaj pkdaUjdpUk + F-[yL + u(r, t) + vt, z] 



(2.24) 
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Here the gradients may be understood either as lattice differences or as the corresponding expressions in momentum 
space. The coefficients are 



d 



= 7% - J2 Q^QjQkQi\UQf^-g^GkiiCl,y^- Q), 
Q ^ 

Ao.j = Y.QrQjQkQi\UQ\h^GkiiQ,^ ■ Q), 



Q 



Ciajpk — 2 X/ QiQiQkQlQm\Uq\'^-^-^Glm{Cl,^ ' Q), 



(2.25) 
(2.26) 
(2.27) 
(2.28) 



where 



_d d_ 



if,,(q = 0) 



is the bare linearized elastic matrix. To leading order, the force is unrenormalized, 



(2.29) 



(2.30) 



Here we have written the expressions in terms of the Fourier transform of the Green's function, G'(q, cj), which is 
defined in the extended zone scheme (i.e. it is periodically repeated in each translated copy of the fundamental 
Brillouin zone). It is 



iuji + K{q) 



9{\q\-A). 



(2.31) 



The 0-function is present because only a partial mode elimination has been performed, so that the slow modes remain 
as dynamical variables in the coarse-grained theory. Note that these expressions diverge in the limit of zero velocity 
and identical periodicities, since in this case all the Q are equivalent to the origin in momentum space. 



2. Random potential 

In the case of the random potential, we adopt the 
approach of disorder pre-averaging the MSR functional, 
thereby working directly with the variance of the random 
force. This is done purely for technical reasons: identical 
results are obtained by coarse-graining for a fixed real- 
ization of disorder and eventually averaging only physical 
quantities. 

We take the simplest model for disorder, in which U{r) 
is Gaussian distributed and statistically translationally 
invariant, with zero mean and second cumulant 



S 



xdidjt[x - x' + u - u' + v(i - i'), z - z']- (2.35) 

As before, we can coarse-grain by integrating out the 
fast (>) modes. The analog of Eq. |2.22| is 



^eff = 5o[u<,u<]-ln(e-^i)g^ 



(2.36) 



In Appendix B, we compute this average to second or- 
der in Si . This again gives an equation of motion of the 
form, Eq. |2.24|, but with 



[t^WC^(r')]e„s. =r(r-r'). 
The force-force correlator is thus 



(2.32) 



A. 



-d.djVir-v'). 



(2.33) 



E 



fyfc,[x + vi,z]r"Gfc;(r,i), (2.37) 



fy7ci[x + vi,z]iG'fci(r,i), (2.38) 



The averaged statistical weight is readily computed: 



1 



[w[u, u]; 

with S = So + Si, and 



(2.34) 



2 ^ Jzt 



iapj 

fyM[x-f vt, z]r"r''Gw(r,i), 



(2.39) 
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-IE 



ijklrn 



[x + vi,z]r"r^G/,„(r,t), (2.40) 



where we have abbreviated F, 



didj • • • f . Also useful 



are the corresponding expressions in momentum space, 



4 - ■ — 



d 



qtqjqkqiT{ci}i—Gki{q.,^- qt), 



(2.41) 



7ii = l5i 



d 



qzqjqkqi^i(Di-g^Gki{ci,^^ ■ qt), (2.42) 



1 



g.gj-gfcgif (q) — ^GH(q,v • qt), 



if ~ d d 

2 y gigjgfe9i9mr(q) — ^Gfc;(q,v-qt)- 



(2.43) 



(2.44) 



A correction 5Tij{Y) is also obtained to the random 
force correlator. Because it is unpardonably ugly, we 
quote it only in Appendix B. It is, however, straightfor- 
ward to show that the renormalized force-force correlator 
cannot be written in terms of a random potential corre- 
lator, i.e. 



r,,(r) EE f,,(r) -f <5r,,(r) ^ -d,d,T{v), 



(2.45) 



for any function r(r). The difference from the equilib- 
rium form can be accounted for by separating the force 
into two components: 



F:{v)^F'r^-{v)+Fr''-{v), 



(2.46) 



where the equilibrium component F^"^' is the gradient of 
a potential, so that 



[Fp-{v)Fp-{v')] _ ^-a,a,r^-(r-r') 



(2.47) 



The other component we denote a (u-independent) static 
force, with the correlator 

[Fr-{r)Fp-ir')] = g,,S{r - r'). (2.48) 

Clearly, since the correction term 5Tij is small, we have 

r<='5-(r) w f(r). (2.49) 

The static force variance is determined, however, entirely 
by the correction. It can be obtained by integrating 



(2.50) 



Substituting in Eq. (8.25) from Appendix B, all but the 
first and fourth terms are total derivatives and hence in- 
tegrate to zero. After a certain amount of manipulation, 
we find 



gigjgfegi'?mgri|r(q)pGfc„i(q, qt- v) 



X [Gi„(q, -qt- v)-Gi„(q, qt- v)] 



(2.51) 



As expected, this expression vanishes for v = 0, i.e. in or- 
der to satisfy the fluctuation-dissipation relation in equi- 
librium, the random force must be the derivative of a 
random potential. 



C. Transformation to Laboratory Frame 

Up to this point, we have worked with conventional 
displacement fields, defined in the crystal frame. This 
means that each particle in the lattice is labeled by its 
equilibrium position r = (x, z), and that its actual trans- 
verse position is given by 



X = x + vt + u(x, z, t). 



(2.52) 



Most measurements in the systems of interest are, how- 
ever, conducted in the laboratory frame. It is therefore 
advantageous to adopt a description based directly in this 
frame. To do so, we define a new field (p{r, t) to be the 
displacement of the particle which is located at position 
r at time t. Formally, this is described by the implicit 
equation 

<?!),(X, z, t) 0,(x + + u(r, <), z, t) = u,{r, t). (2.53) 

Here, because u is defined only at a discrete set of points 
{x}, we have written the previous equation in terms of a 
smoothed continuum field u defined at all space points. 
For many purposes this distinction is insignificant, but it 
will return at one important juncture. 

At this point the need for coarse- grain ing prior to the 
frame transformation is clear. Eq. (2.53) has an unam- 
biguous solution defining (p only when daUi <C 1. In this 
case, we can obtain the transformation rules for gradients 
by simply differentiating. They are 



dtUi = (Sij + dj(j)i) ((9f + V • V) </>j + 



daUi = da(f)i + da(f>jdj(t)i 



dadpUi = dadi3(t)i + 



(2.54) 
(2.55) 
(2.56) 



We are now i n a p osition to transform the equation 
of motion, Eq. ( 2.24 ). The first step is change from u 
field to the u field. This is done by multiplying by a 
smoothing function and summing, using 



u(x, z, t) = a'^' ^ L'(x - x')u(x', z, t), 



(2.57) 



where D{x.) is a delta-like function smoothcd-out on the 
scale of the lattice spacing a, i.e. / , -D(x' — x)/(x') = 
/(x) and D{0) = 1. 

Carrying out this procedure, the gradient and time- 
derivative terms are essentially unchanged, with u — s- u 
to a good approximation. The discreteness of the lattice 
sum is important for the force term, however. It becomes 
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D{x - x')-F/[x' + vt + u(x', z, t), z] 

x' 



Q 



i^"[x + vi + u(x,z,t),z]. (2.58) 



Making the final transformation from u to (p, the force 
term becomes 

i^,[X,z;0,t] = ^e*'5-(^-^*-'^(^^^'*»i^/[X,z]. (2.59) 
Q 

Putting this together with the gradient transforma- 
tions, Eq. ( 2.24 ) becomes 



Iijdt4>j0^i t) — Aiajda4>j + BiapjdaOpYj 

+C^aj(3kda4>jdp(j)k + F,\X, z; 0, t]. (2.60) 

At this point, we can regard X as a dummy variable and 
treat Eq. 2.60 as simply a continuum partial differential 
equation. The transformed gradient coefficients are 



7^J 


= 7ij, 


(2.61) 


4 ■ ■ 




(2.62) 






(2.63) 




= Ciajf3k — lij^ak'^^ + Mak^'j^. 


(2.64) 



III. LINEARIZED THEORY 

In the previous section we derived the p roper form 
of the hydrodynamic equations, Eqs. l.ll| , describing 
the low-energy, long-wavelength properties of a periodic 
medium driven through either a periodic or quenched 
random potential. These equations are distinct from 
their equilibrium counterpart in three respects. First, in 
addition to the well-known convcctive term arising when 
transforming from the crystal to the laboratory frame, 
the coupling of the driven system to the external poten- 
tial yields several other terms linear in the gradients of 
the displacement field; these are collectively described 
by non-zero coefficients Ai^j . Second, the equations con- 
tain nonequilibrium KPZ-type nonlinearities that can be 
thought of as corrections to linear elasticity. Third, in 
the random case there is a nonequilibrium component of 
the static (i.e., (/)-independent) pinning force. This force 
is a genuine nonequilibrium effect, as it cannot be repre- 
sented as the gradient of a potential and vanishes in the 
absence of external drive. 

The remainder of this paper focuses on the random 
case, which is of more immediate experimental relevance. 
Although, as discussed in the Introduction, a general 
analysis of Eqs. 1.11, is beyond the scope of this paper. 



which only the leading (in F) contributions to each term 
are kept. In this approximation the damping = jSij 
is diagonal, and the leading linear gradient term is sim- 
ply the convective derivative, Aiaj —"fSijv". Similarly 
keeping only the bare elastic matrix, we obtain the sim- 
plified set of nonequilibrium elastic-hydrodynamic equa- 
tions for the driven lattice, 



+F,[r,(f},t]+ij,{r,t) 



(3.1) 



where is the usual elastic matrix, rji represents "ther- 
mal" noise (which as we will see in Sec.^ can appear even 
at zero physical temperature) , and the random force Fi is 
given by Eqs. 1.12 and 1.13 and contains all nonlinearities. 



Note that we have kept both the equilibrium and non- 
equilibrium components contained in Fi, since, although 
these are of differing order in F, they cannot be regarded 
as corrections to any zeroth order terms. On e important 
ingredient which has been left out of Eq. 3.1 is the set of 
KPZ non-linearities [Ciajpk terms). While these coeffi- 
cients are probably small, they may very well modify the 
asymptotic long-distance behavior. Based on previous 
work on somewhat simpler (but instructive!) modelsJH3'E3 
we expect that these effects will only increase the distor- 
tions of the moving lattice, and furthermore, that the 
increase in displacements will be most pronounced in the 

longitudinal direction. 

We therefore proceed with the analysis of Eq. 3.1, ex- 
pecting that our results underestimate the roughness in 
the moving lattice, and therefore provide a necessary (but 
not sufhcient) condition for its stability. To proceed, sue 
note that, as pointed out by Giamarchi and LeJDoussalES, 
many of the terms in the random force, Eq. 2.59, are os- 
cillatory in time. Such oscillatory terms average out at 
large sliding velocities, when the "washboard" frequen- 
cies lvq = Q • V are large. Even when they are not large, 
they generate only finite renormalizations of the other 
parameters in the model. We therefore drop this oscilla- 
tory part of F, keeping only the Fourier components of 
the force orthogonal to the mean velocity. 



J2 e'"^-^^- 

Qv=0 



^^F:{r), 



(3.2) 



where we have dropped the explicit time dependence t 
inside Fi, since it has disappeared in this approximation. 
Recall that the (/)-independent force (r) contains both 
equilibrium and nonequilibrium contri butions, w ith the 
corresponding correlators given in Eqs. 



2.46-2.51 



we take this section to discuss the weak-disorder limit, in 



In this section we will further simplify the problem 
and linearize the nonequilibrium hydrodynamic equa- 
tions (3T). We will then discuss the predictions of the 
linearized theory for the decay of translational and tem- 
poral order in the driven lattice. 

We recall that we have used the labels longitudinal 
(l) and transverse {t) to denote the directions along the 
oriented manifolds and transverse to them, respectively. 
The d-dimensional position vector was then written in 
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terms of a transverse coordinate x and a longitudinal 
coordinate vector z as r = (x, z) . We now choose the x 
direction (one of the dt x coordinates) along the direction 
of the driving force (v ~ wx) and denote by r_L the d—1 
directions transverse to the external drive {d = di + dt). 
These d—1 directions can be further broken up into 
components transverse and longitudinal to the oriented 
manifolds, r_L = (y, z), with y a dt — 1-dimensional 
vectors denoting the components of the transverse dis- 
placement that are also perpendicular to the external 
driving force. The di-dimensional transverse coordinate 
space is then described by a 1-dimensional coordinate x 
(not to be confused with x) and y, respectively paral- 
lel and perpendicular to the direction of motion, with 
X = (a;,y). The d-dimensional position vector is written 
as r = (x, z) = (a;,y,z). We will also use the labels par- 
allel (II) and perpendicular (_L) to denote the directions 
parallel and perpendicular to the external drive. For con- 
creteness we specialize to a lattice where the di directions 
along the manifolds that compose the lattice are isotropic 
and the longitudinal elastic properties are described by 
a single elastic constant, denoted by C44. We assume the 
df-dimensional lattice is described by isotropic elastic- 
ity, with two elastic constants, a compressional modulus 
cii and a shear modulus cee- Our model hydrodynamic 
equations for the driven lattice are then given by. 
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+ciidiVt ■ct) + F,[r,cp], 



(3.3) 



where for the purposes of this section we have dropped 
thermal noise, thereby ignoring the suhdominant ther- 
mal fluctuations. These equations describe the impor- 
tant physical case of the lattice of magnetic flux lines in 
a three-dimensional superconductors {di = 1, dt = 2). 
For di — Q and dt = 2 the equations describe the elastic 
properties of driven vortex lattices in superconducting 
films, driven magnetic bubble arrays, or driven Wigner 
crystals. 

The displacement field can be split into components 
parallel and perpendicular to the mean motion, (p ~ 
{<j}xT<t>y)- The transverse displacement <f)y is a {dt — 1)- 
dimensional vector. The pinning force is then seen to be 
independent of (p^'- 

Qv=0 

x[^^;='^(r)+i^r'^(r)], (3.4) 

The independence of the static pinning force on (p^ is a 
consequence of the precise timc-translational invariance 
of the system. This requires that the equation of motion 
be unchanged upon transforming t ^ t+r, (f)^ if)x+vT. 
As argued above, all explicitly time dependent terms are 
irrelevant at low frequencies, thereb y im plying indepen- 
dence of (px in the same limit. Eq. j3.3| is thus linear in 
the longitudinal displacement (px , which can therefore be 
treated exactly. 



The o nly r emaining nonlinearity in the equation of mo- 
tion, Eq.3.3 is in (py, entering through the random force 
Fi[r,<py], which, in this section, we treat perturbatively. 
We stress that the validity of this perturbative calcu- 
lation requires the displacements <py transverse to the 
mean motion to be small, but places no constraints on 
the size of the displacements along the direction of mo- 
tion. The perturbation theory in Fi [r, cpy] thus gives a 
complete description of the behavior of positional corre- 
lations in the x direction out to asymptotic length scales 
for the model of Eq. 3.1. In contrast, the predictions of 



the perturbation theory for the correlations of the trans- 
verse displacement cpy presented below are o nly va lid in 



must 



the Larkin regime; RG methods of Sections 
be used to go beyond the Larkin length scale. 

To lowest order in the random force, we sim ply neglect 
the displacement field <py in Fi[r,cj)y], Eq. (3^). The 
Fourier components of the pinning force are then given 

by 



Qv=0 



(3.5) 



with correlations 

= (27r)'^'<5('^')(q, + qi)(2^)'^'j('^')(qt + q;)Ay, (3.6) 
and 

Qv=0 

= (A + go){Stj - 5ix5jx) + giSixSjx- (3.7) 

The disorder strength A is the variance of the equilibrium 
part of the static pinning force. 



(3.8) 



Qv=0 



The coefficients and gi are determined by the correla- 
tions gij of the nonequilibrium part of the static pinning 
force, according to 



gij — goi^ij — SixSjx) + giSixSjx- 



(3.9) 



They are evaluated in Appendix C, where it is shown 
that for large sliding velocities and short-ranged pinning 
potential (f ^ a, with ^ the range of the pinning poten- 
tial). 



A2 



5o,i 



y(2df~di)/2- 



(3.10) 



The components of the mean square displacement ten- 
sor are given by 



(3.11) 



B,,(r) =2 / / [l-cos(q.r)] 

x[(0,(q,^))(^*(q,^))]^ 
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The correlation functions of the displacement field are 
easily calculated in Fourier space. The linear equation of 
motion for the disorder induced displacement yields 



The asymptotic behavior of the scaling function is given 

- <,(3-dO/2 for 



(3.12) 



where GL{q,uj) and GtCq, i^) are the longitudinal and 
transverse elastic propagators, respectively, 



GL(q,w) = - 
GT(q, w 



1 



i-f{uj - vq^) + (cii + C66)g? + C44g2 ■ 

1 



(3.13) 
(3.14) 



and PMqt) = quqtj and P^Aqt) 



' quqtj are the fa- 



miliar longitudinal and transverse projection operators, 
with qt ~ qt/qt- The correlation functions of the dis- 
placement field are given by 



[KMq,^))? 



27r(5(w) 



91 



qy 



[|(0y(q,c.))fl_ = 27r<5(c.) 



5(0), 



A + 50 



(3.15) 



{ivq^Y 

2 



+0(5i%|GtP 



' [(Cll -t- C66)gt + C44'7z]^ 

(5(0), (3.16) 



where 5(0) = (27r)''+i5(tJ = OW*(q = 0). Because the 
wavevector integral in Eq. |3.11 is dominated by q^ 



{cql/iv), terms containing (g^/g^)|GTp or {ql / ql)\G 
yield less divergent (generally bounded) contributions to 
the fluctuations compared to the terms without this an- 
gular factor. The behavior in real space is displayed here 
for the case = 0, relevant to vortex lattices in thin su- 
perconducting films. The corresponding expressions for 
di = 1, describing flux- line arrays in three-dimensional 
superconductors, are given in Appendix D. For di — 
we obtain, 



Bxx{x,y) - gi 



\y 



3-dt 



■J n 



{dt) ( F|C66 



(3.17) 



where 



d'^'-iu 1 r 



7w(cii -I- cge) 

{dt)(\x\{cil +C66) 



jvy 



(3.18) 



(27r)° 



1 - cos(y • u)e-"" . (3.19) 



by J^q'^'''(0) = constant and 
s > 1. 

Fluctuations in the direction of the driving force are 
dominated by the nonequilibrium part of the random drag 
and by s/iear modes. As discussed earlier, the correspond- 
ing power-law scaling of 0^ holds out to arbitrary length 
scales. Since 



for > (c66/7^^)|2:| : 



(3.20) 



S.. ~ Ixl^^-''')/^ , for <s: {cee/jv)\x\ , (3.21) 



longitudinal density correlations are short-ranged in dt < 
3, with a stretched-exponential decay. Since, as argued 
above, this behavior will persist even once the full nonlin- 
ear form for Fi[r,<py] is taken into account, the validity 
of the elastic model itself is in doubt. In particular, for 
the physical case of dt = 2, the spatial correlations are 
in fact exponentially decaying, and it is natural to ex- 
pect that this drives the unbinding of dislocations with 
Burger's vectors parallel to the drive. If this is the case, 
both translational and temporal order for Q • v 7^ are 
short range in the driven steady state. 

In contrast, both the equilibrium and nonequilibrium 
parts of the static pinning force contribute to fluctua- 
tions in the transverse direction, which are controlled by 
the compressional modes of the system. For large sliding 
velocity or weak disorder and intermediate length scales, 
the equilibrium part of the static pinning force will domi- 
nate as go ^ (A/w)'^ ^ A in this limit. The perturbation 
theory employed here breaks down, however, at lengths 
larger than the Larkin lengths, Kj! and defined by 
B±±{R'^, R^) ~ where ^ is the range of the pinning 
potential. The Larkin volume is defined for deformations 
of the lattice in the transverse direction (d 



is anisotropic, with 



'C^wcii\i/(3-dt) 



^..(ifll) 



5) and it 

(3.22) 
(3.23) 



It is elongated in the longitudinal direction at large slid- 
ing velocities. Longitudinal displacements ((px) grow 
without bound on all length scales in dt < 3 and a Larkin 
domain cannot be defined in this case. As usual, the 
nonlinearities in the static pinning force must be incor- 
porated non-perturbatively to describe the asymptotic 
decay of the correlations beyond the Larkin length scales. 
We will consider this problem in more detail in the fol- 
lowing sections. A natural suspicion is, however, that 
at asymptotically large distances at a finite temperature 
in two dimensions, the nonequilibrium random drag, go, 
will dominate, destroying also the transverse periodicity. 

We stress that longitudinal fluctuations - that grow 
without bound - are induced by shear modes of the lat- 
tice. This is consistent with the intuition that shear 
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modes (rather than compressional ones) play the dom- 
inant role in melting a crystal. In contrast, trans- 
verse fluctuations are controlled by compressional modes. 
Compressional modes are present in both solids and liq- 
uids and are generally not expected to generate the un- 
bounded strains needed to yield dislocation unbinding. 




Fig.4:Schematic illustration of a driven line smectic, i.e. 
smectic in three dimensions, with dt = 2, d/ = 1. 



IV. DRIVEN SMECTIC DYNAMICS 

It was shown in the previous section that disorder- 
induced fluctuations parallel to the direction of the mean 
motion grow algebraically in dt < 3. This algebraic 
growth of fluctuations yields short-range positional cor- 
relations along the direction of the external drive, and it 
is our belief that this implies the breakdown of the elas- 
tic description along the direction of motion (x) . While 
we have not considered explicitly the role of dislocations 
here, analysis _a£_similar equations of motion including 
non-linearitie£3cB suggests that this should occur (for 
dt < 2) via the unbinding of dislocations with Burgers 
vectors parallel to the driving forceE3. This mechanism 
should convert the longitudinal translation and tempo- 
ral correlations to the exponential (or stretched expo- 
nential) form typical of a liquid. This is also in agree- 
ment with real space images of driven two-dimensional 
vortex lattices Idt = di — 0) obtained via numeri- 
cal simulations^ Therefore for dt = 2 and di — the 
driven lattice can only retain periodicity at reciprocal 
lattice vectors transverse to the direction of motion. At 
best, therefore, it consists of a stack of one-dimensional 
liquid-like channels, sliding parallel to the direction of 
motionll3 and has the spatial symmetry of a smectic liq- 



uid crystal. It is, however, important to stress that in 
real, equilibrium smectic liquid crystals the underlying 
rotational symmetry (which is broken spontaneously) en- 
forces soft Laplacian in-layer elasticityt^ In contrast, the 
rotational symmetry of the nonequilibrium "driven smec- 
tic" discussed here is explicitly broken by the external 
drive. These systems are therefore characterized by con- 
ventional gradient elasticity. 

For the case dt = 2 and di = 1, corresponding to mag- 
netic flux-line arrays in three-dimensional superconduc- 
tors, the driven crystal is potentially stable as transla- 
tional correlations decay logarithmically also in the di- 
rection of the applied force. We recall, however, that 
the results obtained here by dropping the KPZ terms are 
expected to strongly underestimate the growth of longi- 
tudinal fluctuations of a moving crystal. Since in d = 3 
the transverse line smectic state is equally stable it is 
possible to have a phase transition between the trans- 
verse smectic (expected at intermediate velocities) and a 
moving glass (stabilized at high velocities). The nonequi- 
librium smectic state (illustrated in Fig. 4 then consists 
of liquid-like sheets of flux lines lying in planes parallel to 
the zx plane. These sheets are periodically spaced in the 
y direction, normal to the external drive. Within each 
sheet, however, there is no positional order of the flux 
lines and the correlations are liquid-like. 

In a smectic liquid crystal the density fleld p becomes 
an independent variable, as it is no longer slaved to the 
displacement field. The density is a conserved quantity 
and therefore a density fluctuation relaxes at a rate that 
vanishes in the long-wavelength limit, i.e., it is a hydrody- 
namic variable. In addition, the system still has broken 
translational symmetry in the direction perpendicular to 
that of the mean motion. Both the Goldstone modes of 
the broken translation symmetry (the displacement (j}y) 
and the conserved density fleld {p) must be retained in 
a hydrodynamic description. In the remainder of this 
section we consider the case dt = 2 and d; = of a two- 
dimensional lattice driven over a disordered substrate. 
The hydrodynamic variables of the driven smectic are 
then a conserved density field p and the one-dimensional 
displacement vector 0^(x, z) (/)(x, z), describing dis- 
placements of the layers in the directions normal to the 
external drive and to the layers themselves. As we are 
only interested in overdamped systems here, the momen- 
tum and energy are not conserved and therefore need not 
be included explicitly in our hydrodynamic description. 

The continuum hydrodynamic free energy for the over- 
damped smectic is given by 




+i^i(a,0)2 + 2if2(9,<^) — |, (4.1) 

where Sp = p — po, with pq the equilibrium density. Here 
cl is the smectic bulk modulus, c\^ the in-layer com- 
pressibility and Ki the layer bending stiffness. The cou- 



15 



pling constant K2 also has dimensions of an elastic con- 
stant. The hydrodynamic equations of the driven smec- 
tic contain additional nonequilibrium terms, as compared 
to their equilibrium counterpart. The nonequilibrium 
terms can be constructed by preserving the invariance 
under inversions about the direction of the external drive 
{y —f —y, (j) —(j)) and the fact that physical properties 
are invariant under the "phase shift" (/) — > (/< + a. 

Density conservation requires that p satisfies a conti- 
nuity equation, 



dtSp + V ■j = 0, 



(4.2) 



where j is the number current density. The equation for 
the layer displacement is given by 



(4.3) 



with Fq a kinetic coefficient and Fy the y-component 
of the non-oscillating part of the pinning force given in 
Eq. 3.2. We neglect here the oscillatory contributions to 
the pinning force that only give small corrections to the 
non-oscillatory part. We also neglect other contributions 
to the pinning force that couple to the density. These 
would be important for a full systematic RG treatment, 
which we do not attempt here. A discussion of more 
general equations for the driven smectic that incorporate 
these terms is left for future work. 

To close the equations, a constitutive relation for the 
current flux j is needed. This is given by 

jx ^ vdp + viSp + poV2dy(l) - poTid^—^, (4.4) 

dp 



jy = PQV3dx4' - P0^2dy 



Sp 



(4.5) 



The first term on the right hand side of Eq. (4.4) arises 
from the transformation to the laboratory frame. The 
terms proportional to the coefficients wi, V2 and V3 are 
nonequilibrium terms that can be generated upon coarse- 
graining, by the method described in Section II for the 
the driven lattice. Other nonequilibrium terms that yield 
contributions that are higher order in the gradients, and 
therefore subdominant, have been neglected here. By 
insert ing the c onst itutive equation for the current in 
Eqs. (0|) and (0|), one obtains, 



(dt + iA)Sp = ~poi2d..dy(f> + [Didl + D2dl]5p, (4.6) 



and 



(dt + i3dx)<l> = D^dy^ + [D^d'^ + Did']<P + — [x, 0]. 

Po " Po 

(4.7) 

Additional "velocities" vi, V2, vs have been defined as 
vi — V + vi, (4.8a) 



V2 = V2 + W3, 
V3 ^ V ~V3 



(4.8b) 
(4.8c) 



The coefficients Di have dimensions of diffusion constants 
and are given by 



Di = TiCi/po, 
D2 = T2Cl/pq, 
D3 = ToKi/pq, 

= (Toc^i - T2K2)/po, 
D5 = (ro/v2 - r2Ci)/po. 



(4.9a) 
(4.9b) 
(4.9c) 
(4.9d) 
(4.9e) 



By solving the hydrodynamic equations in the long 
wavelength limit, one finds that the decay of density 
and displacement fluctuations is governed by two diffu- 
sive modes with eigenfrequencies. 



= V3qx - I 
p - ViQx + i 



D3,1 + (d,-J^),1 

\ Vi — V3/ 



Ul 



Diqi 



D2 



Vi - V3 



V2D5 



(4.10) 
(4.11) 



> and D2 



For stability we must have D4 vi-va 

The diffusion constants D^ and D3 are 



V2D5 



V1—V3 



> 0. 



clearly positive defined. The first mode describes long- 
wavelength deformations of the layers and governs the de- 
cay of displacement fluctuations. The second mode cor- 
responds to the permeation mode of the driven "smectic" 
and describes the transport of mass across the layers that 
can occur in these systems without destroying the layer 
periodicity. It is associated with density fluctuations and 
it has important physical consequences for the response 
of the driven smectic to an additional small driving force 
applied normal to the layers (see Sec. VII). 

An important physical quantity that can be measured 
in both simulations and experiments is the structure fac- 
tor of the driven periodic medium. As the driven smec- 
tic has broken translational symmetry in the y direction, 
normal to the layers and to the external drive, the density 
fleld can be written as 



p(x) = PL (x) + pQ^ (x)e' 

Qy 



JyV 



(4.12) 



where Qy = n2Tr/a are the reciprocal lattice vectors cor- 
responding to a layer spacing a, with n an integer, pg 
are the corresponding Fourier components of the density, 
and pl(x) is the smooth (liquid- like) part of the density 
field. The smectic structure factor is then given by 

S{ci±) = SLiq±) 

+ f e-*«--e-*(««-««)^(pQ„(x)p_Q„(0)), (4.13) 

Qy 

with S'L(q±) = {\pl{<\i_)\'^) ■ The smectic structure func- 
tion should therefore consists of a broad liquid-like back- 
ground Sl^SIi.), with superimposed peaks at the recip- 
rocal lattice vectors Qy, normal to the direction of mean 
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motion (see Fig. 2(b)). As discussed in the Introduction, 
in a Gaussian theory the correlator of the order parame- 
ters pQy (x) can be written in terms of the mean square 
displacement, according to 



p2g-Q^([0(x)-0(O)]^>^ 



(4.14) 



The disorder-induced transverse mean square displace- 
ment is easily calculated from Eqs. 4.6 and 4.7 by treat- 
ing the the random force ^^/[x, (j)] as a perturbation, in an 
analysis similar to that of Sec. [11. Conceptually simple 



but algebraically tedious calculations show that including 
the coupling to the density does not change the resulting 
decay of the correlation function. For s implicity we there- 
fore neglect this coupling in Eq. (4.7) (i.e., let 5p — 0) 
and obtain 



[m^)-m?)]. 



(A 



x\Da 



War 



(4.15) 



The scaling function (s) is identical to that ob- 
tained in Eq. (3.17) and it has the asymptotic behavior 



= constant and T^ '{s) ~ for s > 1. The 
perturbation theory therefore predicts that the smectic 
Bragg peaks at qy — Qy decay exponentially with the 
system size. In other words, disorder would destroy the 
transverse periodicity of the smectic. In fact, we will see 
in the following sections that this result continues to hold 
non-perturbatively. The power-law scaling of transverse 
Bragg peaks obtaincji-in simulations of two-dimensional 
driven vortex latticeaHj is therefore most likely an artifact 
of small systems and weak disorder, and would crossover 
to a disordered form at longer distances. We will return 
to this point later in Sec. VII. 



(77(r,<)ry(r',0) = 27T5(r-r')5(^-^')■ 



(5.2) 



The random force F is characterized by the correlator 
[F(</., r)F(0', r')]_. = A(0 - 0')<5(r - r')- (5.3) 

The function A{(j)) is periodic with the smectic lattice 
spacing, which we take to be a = 27r. It includes both the 
equilibrium and non-equilibrium components. The latter 
can be viewed as simply an overall constant contribution 
to A(0) . In addition to discounting density fluctuations, 
Eq. 5.1 also neglects an allowed KPZ non-linear term of 
the form 



"KPZ 



(5.4) 



which should be added to the right hand side of the equa- 
tion of motion. Simple power-counting shows, however, 
that, in contrast to a driven lattice this term is strongly 
irrelevant in a tran sver se smectic (see below) . 

To analyze Eq. (|5.l| ), we use the method of MSR to 
transform the stochastic equation of motion into a field 
theory, similar to what was already done in the Sec. [IB 
to perform the single-step coarse-graining. In this case, 
the MSR "partition function" is 



Z 



[dcf)] [d(j)\ e 



(5.5) 



where S = Sq + Si, with 



(5.6) 



and the interaction term 



V. RENORMALIZATION GROUP FOR "TOY" 
SMECTIC 

A. Model and MSR formulation 

In this section, we will consider a simplified model for 
the smectic phase, in which the hydrodynamic fluctua- 
tions of the conserved density are neglected. Th is " toy" 
smectic is then modeled simply by dropping Eq. 4^ and 
setting 5p = const., leaving the single equation of motion 



l{dt + vd,)cp ^ Kud^cj) + A'iVi0 + F{(^, r) + 



(5.1) 



Here we have pulled 
if II = Dspo/To, 



out a factor of 7 = Pq/^Oj l^t 
= Dipo/To, and F± F. For 
simplicity we assume longitudinal and transverse elastic- 
ity are governed by the same elastic constant Arj_. We 
have also added the random time-dependent "thermal" 
noise ry, satisfying 



Si 



(5.7) 



By construction, Z = \; non-trivial correlation and re- 
sponse functions arc obtained, however, by inserting ap- 
propriate (j) and (p operators into the functional integrand. 



B. Power counting 

Let us first consider under what conditions the random 
force is a relevant perturbation in the sense of the renor- 
malization group (RG) using simple power-counting. To 
do so, we rescale the coordinates and fields by a scale 
factor 6 > 1: 



xj_ &xj_, 
X b'^x, 
t -> b^t, 

4) b^4). 



(5.8) 
(5.9) 
(5.10) 
(5.11) 
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Note that, anticipating a periodic fixed point, we will 
not rescale (f). To fix the exponents ^ and x, we choose to 
keep the terms and K± fixed in the quadratic action 
Sq. This implies d~3 + z + ( + x = and d~l + z + x — 0, 
or 



c 

X 



1 



(5.12) 
(5.13) 



Note that because was kept invariant, the temperature 
T necessarily rescales 



(5.14) 



Clearly this is a somewhat artificial choice for A = 0, 
since the theory retains no memory of the periodicity of 
(j). It is rather natural, however, for A ^ 0, and will 
be returned to below. The exponent z is more subtle. 
Naively, it should be determined by the condition that 
7 be invariant under rescaling. Neglecting the random 
force, we then obtain 



= 2. 



(5.15) 



We will see that even for small non-zero A, this is actu- 
ally very far from correct. 

Fortunately, the rescaling of A(0) and C is in fact in- 
dependent of z. The KPZ non-linearity is strongly irrel- 
evant: 



C-^ 0, 
while the disorder correlator obeys 
A(0) ^ 63-'^A(0). 



(5.16) 



(5.17) 



We therefore see that for d > 3, the random force is irrel- 
evant, and at long length scales has a negligible effect on 
the moving smectic. For d < 3, however, power-counting 
is insufficient to determine the fate of the system. 



C. Zero temperature RG 

To proceed, we will perform a systematic RG, working 
perturbatively in Si, expanding it from the exponential 
and integrating out "fast" modes with large transverse 
momenta. This is simplest in the scheme with neither 
a frequency lo nor longitudinal momentum cut-off. 
The fields are first decomposed into their slow (<) com- 
ponents with q±^ < A/6 and fast (>) components with 
A/6 < g_L < A, where A is a hard transverse momentum- 
space cut-off. That is. 



Then the partition function is 



(5.18) 



Z= /[d(^<][d<^<]e- 



-5o[0<,0<] /„-Si[0<-i-<;'>,0<+0>] 



(5.19) 



where the brackets denote an average over the fast fields 
with respect to the quadratic action S'o[0>,0>]. Rela- 
befing the surviving slow fields <^< ^ 0, 0< — > 4>, the 
renormalization of the effective action due to the mode 
elimination is given by 



-Sett. 



-Si[<p+4>>,(i>+<t>>] 



A cumulant expansion gives 
(e-^0>=exp|-(^: 



'i/> 



(5.20) 



(5.21) 



where the subscript c indicate the cumulant (connected) 
correlator. From Eq. 5.14, it is natural to suspect that 
much of the interesting physics is dominated by small 
(renormalized) temperatures. With this in mind, we first 
consider the RG in the extreme T = case. Then the 
only non-vanishing expectation values correspond to re- 
sponse functions: 



G>(r,i) = (0(r' + r,i' + t)<^(r',i')>: 



(5.22) 



A/6 



dqx 
2ti 



duj 



27r i^{uj - vqx) + K\\ql + K±qj_ ' 



It is instructive to perform some of the integrations above 
and bring out the dependence on x and t. We are primar- 
ily interested in the limit K^^ 0, since i^y is obviously 
less relevant than "fv, which has one less x derivative. In 
that limit these integrals are trivial, and we obtain 



G>(r,i) = 



A/b 



d^ 

(27r)'i-i 



7 



^-iK^/'y)Atg(^t)S{x-vt). 

(5.23) 



From this we see that the response function is causal and 
represents unidirectional propagation along the positive 
X axis. Including a non-zero i^y simply spreads out the 
5-function over a distance Sx cx y/K^^t. The irrelevance 
of i^ii is indicated by the smallness of this width relative 
to the distance vt propagated in the large t limit. 
Working first to 0(A), consider the term 

0rt[A'(O) + A"(0)(i - i')at</'rt]G>(0, t - t'). (5.24) 



By symmetry, A'(0) = 0, and we need to evaluate 
rG>(0,r) = 



^ i-f{uj - vq^) + K\\ql + K^q\ 
dl 



27/f||A^-iGd_i 



(^2^2 + 4^ Jf_^A2)3/2 



-'d-l 



1 + 2k vK± 



dl, 



(5.25) 
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where we defined 



K = 



(5.26) 



quencc of taking an ultra- local (i.e. 5-function correlated) 
random force. For such a force, it is straightforward to 
show that the static response function 



chosen the infinitesimal rescaling factor & = e'*', and de- 
fined Cd = Sd/{2ny' in terms of the surface area of a 
d-dimensional sphere Sd = 2tt'^^'^ /T{d/2). Inserting this 
result above gives a renormalization of the friction drag 
coefficient (the inverse mobility): 

^ -^^tgtiA"(0). (5.27) 

Note that we have obtained no renormalization of the 
spatial gradient terms. This is actually a general conse- 



G{q,uj = 0) = 



-i-yvQa; + K«ql + K±q1 



(5.28) 



is unrenormalized, i.e. exact even when F is included in 

the equation of motion. Thus 7^, K^\, and K± suffer no 
diagrammatic corrections at any order. 

The next step is to examine the renormalization of A. 
The first corrections (for T = 0) arise at O(A^). To de- 
termine these, we must compute the next term in the 
cumulant expansion. This is 



6S2 



^>riti'Arit;A((^riti - (/)rit', )^r2t2'^r2t^ ^(^r^t^ - (/jr^t^) 



(5.29) 



At this point we are aided by a simplifying feature that actually makes this part of the calculation easier than 
the equilibrium one. Since wo have a propagating mode along the x axis, in the limit of zero longitudinal damping, 
= 0, the response function vanishes unless a; > 0. This means that two response fields at different x points 
cannot be contracted, since this would give the product G(x — x', ■)G{x' — x, •), which vanishes for any x. Instead, a 
non-vanishing contribution obtains only when both ^ fields are taken from the same term. This gives 



^^2 



1 



r2t2t' 



X [G> (ri -r2,ti-t2)G> (ri - ra, - 4) - G> (ri - rs , - t2)G> (ri - rs, t[ - 4)] . 



(5.30) 



Each term contains two response functions, which constrain their arguments to be small. To leading (zeroth) 



order in gradients of </>, we may approximate (/>r2t2 
4>r2t'2 ~ 0rit'j- This leads to the simpler formula 



and 



'r2t' 



in the first term, while in the second 



55s 



= / 4>riti4>rit[^" {(t>riU " ^rit'J [A(0) - A{(f)r,ti " (t>T^tO] ' 



where the same integral / occurs in both terms. In the limit i^n 0, it becomes 

I-- 



f G>(q,a; = 0)(?>(-q,a; = 0) pa ^'^J'^ dl. 
qj- Jqx 



2-yvK± 



(5.31) 



(5.32) 



This gives the second order contribution to the renor- 
malization of A(^): 



diA{cj>)\ 



^''-'A"{ct>) [A(<^)-A(0)]. 



(5.33) 



At T = 0, these arc all the necessary mode-elimination 
contributions. Combining these results with the scale- 
changes, we thus arrive at the zero temperature RG equa- 
tions 



da = 

diK\\=-2K\\, 
diK = —2k, 



K'^-^Cd 



1 + 2k ^vK± 



-A"(0) 



7, (5.35) 

(5.36) 
(5.37) 



where the flow for the dimcnsionless parameter k is exact 
and independent of the choice of the rescaling exponents 
z, C and X- 



/. Behavior in d = 3 — e dimensions 



diA{<t>) = (3 - d)A{ct>) 



.^-^A"{<l>) [A{<t>) - A(0)] , (5.34) 



As a first attempt at an analysis of these results, we 
consider the problem in d = 3 — e dimensions. Defining 
the dimensionless random-force correlator 
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m 



-A(0) 



(5.38) 



the RG flow equation for this dimensionless random force 
variance is given by 



diTi^) = (3 - dm) 



^r"(0) [r(^) 



r(0)]. (5.39) 



For e ^ 1, we thus expect to fin d a fixed point with 
r = 0(e). However, evaluating Eq. 5.39 at = demon- 
strates that (in contrast to equihbrium case, where there 
was an additional stabilizing — r'(0)^/2 term) this is only 
possible if r( 0)* = 0. Like other functional RG equations 
at T = 0, Eq. 5.39 leads to a non-analytic force-force cor- 



relator. This can be seen directly by differentiating twice 
above and evaluating at the origin: 

diVio) = (3 - rf)r"(o) - i [r"(o)]' . (5.40) 

Since T"{Q;l = 0) < 0, this equation leads directly to 
a divergence. In fact, all "flxed points" (see below for 
an explanation of the quotation marks here) have a slope 
discontinuity at the origin. A little analysis demonstrates 
that the large I behavior of T is 



mi) 



IT — 27rn) 



ro(0, (5.41) 



where Tq is the zero Fourier component of T((j)) and sat- 
isfies 



(5.42) 



Note that although Tq = — 27r^e/3 provides a formal fixed 
point solution, it is unphysical: for this value, r(0) = 0, 
which would imply the (positive semidcfinite) variance 
of F{<j),Y) vanishes, and hence that F is uniformly zero, 
clearly in contradiction with r((/)) ^ 0. For any phys- 
ical situation Fo will be larger. For instance, equilib- 
rium initial conditions require that the force arise as a 
derivative of a random potential, and hence that the in- 
tegral over A((/i) vanish. For the solution above, this 
gives Fo = > — 27r^e/3. Indeed, Fq can be identified 
as the variance of the static for ce ca lculated using the 



single-step coarse graining in Sec, II B 



For physical situations, it is clear from Eq. (5.42) that 
Fq is a strongly relevant variable which will fiow off under 



C>(r,t) 



[>T+r' ,t+t' iPr' t' / > — 



the RG. Luckily, as pointed out 1^ Narayan and Fisher 
in their study of CDW depinningE3, this does not really 
present a problem. It can be easily shown that the static 
random force does not effect the dynamics, by the same 
argument used above to demonstrate the exact static re- 
sponse function (i.e. shifting it away and showing that 
the distribution of the shifted force is unmodified). This 
is also essentially the same argument used in the Cardy- 
Ostlund problem, in which there is also a runaway ran- 
dom force (which is, in that case however, the gradient 
of a potential). 

The non-analyticity of F has important consequences 
for the dynamics. Indeed inspection of Eq. 5.35 shows 
that it is problematic: the quantity A"(0) is (minus) in- 
finity at the fixed point (and indeed becomes infinite at 
a finite length scale) . Physically, the non-analyticity of F 
is related to the existence of multiple metastable minima 
in the effective potential on the scale I and a concurrent 
"sharpness" of this potential. At zero temperature, this 
sharpness leads to trapping of the phase </> and hence 
to a breakdown of the simple assumption of analyticijty, 
in the coarse-grained dynamics. A detailed analysisE3 
shows that the divergence of A"(0) can be interpreted as 
the signature of a multi- valued force reminiscent of static 
friction, and corresponds to the existence of £i_cion-zero 
transverse critical force to "depin" the smectidEl. 



D. Finite temperature RG 

A simpler means of controlling these singularities is to 
include non-zero thermal fiuctuations, which act to lo- 
cally average the effective potential and effect thermally 
activated motion between different metastable states on 
long time scales. In hopes of obtaining a workable dy- 
namics, we are thus led to consider the effects of a non- 
zero temperature. At the same time we must also ask 
if temperature is a relevant or irrelevant perturbation 
around the T = fixed point considered here. Naively, 
from Eq. 5.14, we would expect T be strongly irrele- 
vant for d near three. In equilibrium, the corresponding 
(power-counting there gives 2 — d) result is exact; again 
due to Galilean invariance and thepFDT, there are no 
diagrammatic renormalizations to Te3. 

Because of the lack of FDT, we will, surprisingly, be led 
to a completely different conclusion! Consider the mode 
elimination at a non-zero temperature. The correlation 
function. 



2jTe 



-ici-r-\-iujt 



72(w - vq,,)^ + {K^l + K^ql)^ 



(5.43) 



is now non-zero. This now makes it possible to contract two (j) fields when renormalizing the MSR functional. It is 
sufficient to work only to linear order in A (or F). Following the same method as earlier, we find 



(T) 



1 



lrt$rt'A"icj)rt " M [C> (0, 0) - C>{0,t-t') 



(5.44) 



This gives two terms. The first contributes a renormalization of A. That correlation function gives 
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c>(o,o) 



2T7 



(5.45) 



Note that the integral is singular in the limit. 

Physically, this is because some damping is needed to 
control the thermal fluctuations excited in the propagat- 
ing longitudinal mode. Incorporating this piece in the 
RG equation for T gives 



dim 



where 



(5.46) 



T = A''-^Cd-i 



T 



(5.47) 



The second term in Eq. 5.44 gives a correction to the 
temperature. This is simple to see, since in this term 
|i — i'l is kept small, so that 0rt' ~ (f'n (and likewise for 
(f)). This gives 



SS 



i £ (0rt)' A"(0) C>(0, t-t'). (5.48) 



This integral is finite in the limit i^n 0, giving 

»> 



T 

K I V 



K'^-^Cd-idl. 



(5.49) 



This is a renormalization of ^T. Taking i nto a ccount the 
renormalization of 7 obtained earlier (Eq. 5.35), one finds 



diT = 



I 1 n 

l-d- [ 

2 1 + 2k 



r"(o) 



T. 



(5.50) 



Using the definition of the dimensionless temperature, 
this becomes 



2-d- 



1 



1 + 2k 



r"(o) 



T. 



(5.51) 



1. d — 3 ~ e redux 



Eqs. 5.46| and 5.51 complete the modified set of RG 
fiows at non-zero temperature. Let us focus again on the 
behavior for d = 3 — e, discussed above for T — 0. As sus- 



pected, the presence of the diffusion-like term in Eq. 5.46 
indeed acts to smooth out the cusp in T{(j>) (a simple 
heuristic argument for this rounding is given in Ref. 34 ) . 
To see thi s, con sider the "adiabatic" approximation in 
which Eq. 5.46 is solved for fixed non-zero T, ignoring 



corrections to diT arising from the scale dependence of 
T. To carry this out, we search for a solution 



r(0,o-r(o,O-(3-rf)f((/.,r(O), 



(5.52) 
: gives 

the flow equation for r(0, Z), which in turn implies an 
equation for T((j)) in the "adiabatic" approximation. 



where f(0, 1) = 0. Evaluating Eq.^.46 



at 



diTiO,l) 



(3-d)r(o) + rr"(o), 



m 



Mf"(0) + if"(0)f(0). 



(5.53) 
(5.54) 



where /i = T/e <C 1. Multiplying Eq. 5.54 through by 
r'(0), allows one to perform one integral, and thereby to 
solve for r(0) exactly in an implicit form. One finds 



dT 



In/i 



ln(l + r/2^) - f 



-1/2 



2(/), (5.55) 



for < (/) < 2tt. For - 2nn\ > ^^/\\nJl\ (with 
integer n), this reduces to the zero temperature solu- 
tion in Eq. 5.41. It contains, however, a boundary layer 
near (j) = 2Tm. Inside this boundary layer T{(f)) remains 
smooth, and a simple computation finds 



^ ' Tln(e/r) 



(5.56) 



Putting this result back into Eq. ^.51 gives, to leading 
order in e. 



diT 



-T 



21n(e/r) 



(5.57) 



The new term leads to a fixed point at a non-zero tem- 
perature given by 



T 



2|lne| 



(5.58) 



Note that in equilibrium, i.e., in v — > limit, k 00, 
7t/(l -I- 2k) — > 1/2 , and two terms appearing inside the 
parentheses in Eq. |5.5l| — exactly cancel due to the FDT, 
preserving the power-counting result and leaving T ir- 
relevant. Here, for any finite v, they do not, and the 
noise renormalization (first term) overwhelms the dy- 
namic renormalization (second term) due to strong ir- 
relevance of K (or equivalently i^||), thereby destabilizing 
the zero-temperature fixed point! The asymptotic flow 
of T is then given by 



2-d 



-r"(0) 



T. 



(5.59) 
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E. Behavior in three dimensions 



C = i-d + T- 



(6.1) 



In three dimensions, we expect that T flows to zero, 
and the above analysis is invahd. To analyze this case, 
con sider first the behavior of V{(f)) at T = 0. From 
Eq. 5.40 , we see that even for d = 3, r"(0) diverges at a 
finite scale, and the non-analyticity remains. Indeed, for 
T = 0, a solution (presumably the asymptotic attractor 
for an arbitrary initial condition) of Eq. 5.39| is 



b-l-T 



0) 



a{l) min < (0 — tt 

n 

+ro(0, 



2iinf 



where 



dia(l) = 



27r^ 



TT 

T 

(5.60) 



(5.61) 
(5.62) 



Thus the amplitude of the "cuspy" part of the disorder 
correlator decays slowly to zero. 



a{l) = 



a(0) 



l + la{Q)' 



(5.63) 



Employing again the adiabatic approximation for small 
T 7^ 0, one finds 



r"(0,;;d = 3) = = 



i{l)f 



T{l)Ha{l)/T{l)) 



(5.64) 



Putting this back into Eq. ^.51 , the "fixed point" for 
e 7^ now drifts slowly down towards zero. The correct 
asymptotic behavior is obtained simply by substituting 
e ail) in Eq. p3.58 , giving 



2|lna| 2P\hil\ 



(5.65) 



for large I. Note that this is faster than the usual case of a 
marginally irrelevant operator (~ 1/0? but much slower 
than the naive power-counting result (~ e~'). 



has a discrete spectrum (defined in the space of 2tt- 
periodic functions). 



£ cos n(f> — (3 — d — Tri^ ) cos nq 



(6.2) 



For T > 3 — c?, therefore, all the harmonics but the 
constant {n — 0) are irrelevant, and A(0) Aq. For 
(3 — (i)/4 < T < 3 — d, we can study the onset of 
non-trivial random-force correlations by truncating the 
Fourier expansion of A at n = 1: 



A((/)) = Ao + Aicos( 



(6.3) 



Eq. |6.3| displays one distinct advantage over the full func- 
tional RG treatment in the previous section: A{<j)) in this 
limit is manifestly analytic, so perturbation theory is uni- 
formly valid. By contrast, the treatment for d ~ 3, while 
physically reaso nabl e, is con sider ably less controlled. In- 
] or Eq. KM, T"{Oy 



deed, from Eq. 



=-2(d-2)is 

0(1) in this limit, and although r{(f>) is 0(e), it is un- 
clear how the singular derivatives might enter into higher 
order corrections. 

To implement the complimentary, strictly controlled 
app roach for d = 2 + e, we simply insert Eq. O 
Eq. 5.46, dropping harmonics with n > 2. This gives 



(6.4) 
(6.5) 
(6.6) 



where Ao,i = Ao^iCd^A'^^^ / {K±^\v\). We pause to 
point out that Eq. |6.3| is equivalent to the physical ap- 
proximation of the equation of motion, 

7(9t + v5,)0(r, t) = {K\\dl + K^Vl)(b{r, t) + Fo(r) 
+ Fi(r)cos[y-</)(r,i)] +77(r,i) , 

(6.7) 



VI. DYNAMICS OF 2+e-DIMENSIONAL 
SMECTIC 

In this section, we study the toy model for the driven 
smectic near two dimensions. Based on an ex trapolation 
of the results of the previous section (Eq. 5.58), we expect 



with 



that in this case the governing stable fixed point should 
occur at a renormalized temperature of order on e. In 
this regard d = 2 plays a special role, since (see Eq. 5.59 ) 
dimcnsionless temperature T becomes marginal. 

For temperatures of 0(1), the character of the func- 
tional RG (Eq. ^.46 ) is substantially changed. This is 
because the linear operator 



2Ao(5('')(r-r'), 
2Ai5('')(r-r'). 



(6.8) 
(6.9) 



For c omplete ness, we rederive Eqs. 3.4- |6.6| di rect ly from 
Eqs. 6/7-BJ in Appendix E. Note that Eq. |6^ i mplies 
that the random ^-independent drag correlator Aq is al- 
ways generated by the disorder and is relevant, even if 
it is not present in the "bare" equation of motion. Al- 
though Ag runs off to infinity, as discussed in the previous 
subsection, its effects can luckily be taken into account 
exactly, by a simple transformation on the field (j). The 
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remaining flow equations for and T{1) contain three 

fixed points 

Gaussian: T* = 0, = , (6.10) 
Zero Temperature: T* = 0, = 2(3 - d) , (6.11) 
Driven Smectic: T* ^ 5 - 2d, A* = 2{d - 2) . (6.12) 

Of these, only the Driven Smectic is globaUy stable. For- 
tunately, it is also perturbative, and hence controlled, for 
d near 2, and indeed becomes exact in the d — > 2"*" limit. 
Furthermore, the Driven Smectic fixed point that is per- 
turbative near d = 2 appears to smoothly match onto the 
finite disorder, finite temperature fixed point that was ob- 
tained at strong coupling by a functi onal renormalization 
group calculation near d = 3 in Sec. 

are also the 



Equations 6.5-6.6 



log of the Cardy-Ostlund fixed lin 
l-f-l dimensional vortex glass stat' 
rough phase 
dom substrat' 




nequilibrium ana- 



ihat describes the 
and the super- 
crystal surface grown on a ran- 
Because of the lack of fluctuation- 



dissipation theorem for the driven system considered 
here, the Cardy-Ostlund fixed line is destabilized by the 
nontrivial renormalization of temperature. This disorder 
generated thermal renormalization (Ai-dependent term 
in Eq j6.6| ) is reminiscent of the "shaking" temperature, 
discussed by Koshelev and VinokurH Note, however, 
that the "heating" found here is a multiplicative rather 
than an additive effect, as was suggested in Ref. 

In rf = 2 + e dimensions the flow equations have an 
interesting spiral structure around the Driven Smectic 
fixed point, as is illustrated in Fig. 5 




Fig.5.:Renormalization group fiow diagram (for 2 < d — 
2 + 0.3 < 3) in the disorder Ai temperature T plane, 
for a set of initial conditions with Ai(0) = 0.2 and T(0) 
ranging from 0.05 to 0.85, in increments of 0.2 

The structure of the renormalization group flow in the 
physically interesting case of d = 2 is displayed in Fig. 6. 
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Fig. 6: Renormalization group flow diagram (for d — 2) 
in the disorder Ai temperature T plane, for a set of ini- 
tial conditions with Ai(0) — 0.2 and T(0) ranging from 
0.05 to 0.65, in increments of 0.2 

In d = 2 the Driven Smectic flxed point moves down 
to zero disorder and merges into the zero-disorder fixed 
line. Despite the absence of a globally stable finite disor- 
der fixed point in d = 2, we expect nontrivial observable 
effects associated with the interesting RG flows displayed 
in Fig. 6. Qualitatively, a moving lattice at temperature T 
with disorder Ai behaves at long times and length scales 
as a thermal moving smectic with an effective disorder- 
enhanced temperature Toff(r, Ai). 

To determine T^s. , we take advantage of the solvability 
of Eqs. 3.5-3.£ in two dimensions. Dividing Eq. 6.5 by 
Eq. 6.6 gives 



dAi{T)_^l-T Ai 



dT T 
This has the exact solution 

-r-cff 



T 



2 — Toff + f/ToS 



(6.13) 



(6.14) 



where / is given by the initial conditions Ai(^ = 0) = Ai, 

T{1 = 0) = T as 



/ = T[Ai-2 + r] 



(6.15) 



The asymptotic effective temperature is determined from 

off 



Eq. 6.14 by setting A-^ = 0, which gives 



Ta«^l + Jl + T{Ai + T-2). 



(6.16) 



This effective temperature has the following meaning: 
asymptotically, the connected correlation and response 
functions of a moving smectic at temperature T in the 
presence of disorder Ai are given by the noninteracting, 
disorder- free functions with T replaced by Toff(r, Ai). 
Of course the details of the approach to the zero disorder 
flxed line and to this effective temperature will determine 
subdominant g-dependent corrections to the correlation 
and the response functions. 
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VII. NONEQUILIBRIUM RESPONSE AND 
CORRELATION FUNCTIONS 



We now turn our attention to the dynamic response 
and correlation functions of a moving smectic. We will 
consider only mean (disorder-averaged) properties here. 
The mean response function is defined by 



R{r,t) 



' <5(0(r + r^t + tO) 
(5r/(r',i') 



(7.1) 



We similarly define a mean correlation function, 



C(r, t) = [((0(r + r\t + t') - 0(r', t')f)\ . (7.2) 

- \ /J ens. 

Given the results of the renormalization group analy- 
sis of the previous two sections, the long length and time 
asymptotics of these functions can be computed using 
standard matching techniques. For the momentum-shell 
regularization used here, this matching is most directly 
done in momentum and frequency space. Consider first 
the correlation function, which for transverse momentum 
k_L and rescaling factor 1 < 5 < |A/k^| satisfies the re- 
lation 

C(k, oj- {K}) = ¥'^+<+^'C{k^b, Kb^, ujb^; {K{b)}). 

(7.3) 

Here {Xi{b)} denotes the set of running coupling con- 
stants at scale b. The prefactor on the right hand side 
arises from the (conventional) definition of the Fourier 
transform, 



C{r,t) 



ko; 



2[1 



(7.4) 



and the dimensionlessness of (j){r,t). 

To calculate the correlators at long length and time 
scales, we will choose b = A/|k_L| ^ 1 and evaluate the 
right-hand-side of Eq. 7.3. This is simple because the 



rescaled correlator is evaluated at a large rescaled trans- 
verse wavevector |k^|6 = A equal to the uv cut-off, at 
which fluctuations are small, and therefore can be taken 
into account perturbatively, without encountering any 
infa-red divergences. However, the computation of the 
rescaled propagator requires a knowledge of the running 
couplings {Ai(6)}. 

These flows have been studied in Sections IV-V. 
Throughout, unless explicitly indicated otherwise, all 
running couplings without arguments refer to the bare 
couplings, i.e., = Xi{b = 1). It is convenient to choose 
z = = 2 and x — — 1. With this choice, K± and •yv 
are invariant, i.e. 



K^{b)^K^, 
yv{b) — jv. 



(7.5) 
(7.6) 



The remaining parameters behave nontrivially. The 
rescaled tempe ratur e and l ongitu dinal elastic modulus 
follow from Eq.5.47 and Eq. 5.37, respectively: 



2(if^if||)l/2_ 



Kb-^. 



(7.7) 
(7.8) 



Ind — 2, the dimensionless temperature flows to the fixed 
half-line (and is hence non- universal) , while in d = 3, it 
flows logarithmically to zero, i.e. 



(7.9) 



For 2 < d < 3 it flows to a universal value (see Eqs. 5.58 
and 3.12D , which is sadly of only formal interest since 
the 2.5-dimensional smectic is currently experimentally 
inaccessible. 



The mobi lity is more complicated, but using Eq. 5.35 
and Eq. 7.S , we obtain 



7(6) = 7e*('') 



with 



log b 



dl 



l + 2Ke 



1^ c 2/ 



r"(o,o 



(7.10) 



(7.11) 



The integrand in Eq. [7.11 is exponentially suppressed at 
large so that $(6) has a finite limit as 6 — > oo. This 
implies the finite renormalization 



7(6) h- 



(7.12) 



Because, for large v (small 7t) the integral in Eq. 7.11 
is dominated by small I (due to the exponential be- 
havior of k(Z)), the constant a is highly nonuniver- 
sal. Above two dimensions at sufficiently high veloc- 
ities and temperatures (k << 1 and T 3> ag in 3d), 
T flows rapidly to a unique fixed-point, and it becomes 
parameter-independent, with a — 2{d — 2). For d = 2, cr 
depends upon the bare disorder strength even for k << 1, 
due to the semi- infinite fixed line with Toff > 1. More in- 
teresting is the ext reme low-temperature limit, in which 
(see Eqs. 5.56 5.64 ) r"(0) becomes singular. In 3d, this 
limit gives 



Tln(ao/r)' 



(7.13) 



The corresponding low-temperature regime in two di- 
mensions is outside the limits of the controlled RG, but 



we expect a result similar to Eq. 7.13 to hold with, how- 



ever, flQ of order one. Note that in all cases 7^ is strongly 
enhanced as the velocity of the moving smectic is lowered. 

The remaining flow parameter is the force-force corre- 
lator A(0). Again, we quote here the results only directly 
in d = 2, 3. 
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In two dimensions, the correlator is well-described 
asy mpto tically in t he single-harmonic approximation 
(Eq. . From Eq. Ap grows linearly with &, 



Ao(6) b^oo A b, 



(7.14) 



where A > Aq, and for Ai ^ 1 and T > 1/2 can be 
estimated by 



Ao- 



1 



A? 



47rAA'^7t; 2T - 1 



(7.15) 



As can be easily seen from EqJeJ, in d = 2, the first 
harmonic flows asymptotically to zero according to 



Ai(6) b'^ Aifo-C^-'-'-i). 



(7.16) 



To further explore the implications of Eqs. 7.19-7.21 



we now discuss the corresponding expressions in real 
space and time. First consider the long-time limit 
CEA(r) = linit^oo C(r, t), 



C'EAiv) 



2A [1 



(7.22) 



This correlator is analogous to the Edwards-Anderson 
correlator in a spin-glass, and represents a static but 
random conformation of the (j> field that persists at long 
times. Performing the above integration, we find two lim- 
its. For d < 4 and x < K^^/jv, the static roughness scales 
isotropically (similarly to the equilibrium case), with fi- 
nite anisotropy arising due to the difference between K± 
and Ki\ 



In three dimensions, we have ins tead l ogar ithmic flows, 
and simple manipulations of Eqs. |3.60| - ^!6^ give 



, I) fc^oo 2TrjvKj_ 



ao 



1 + Iuq n 



in { (0 — TT — 2?™)^ ~ ""^ } 
(7.17) 



where 



A — 2ttjvK± 



2^2 



-ao 



(7.18) 



We are now in a position to evaluate the right-hand- 
side of Eq. 7^. Setting b = e' = A/|k^| and using the 
above results, we obtain 



where 



I^^^Sf^ + A<5(.) 



d = 2 
rf = 3 
(7.19) 



V{k,iu) = {-iRU-^vh,f + (K\\kl+Ki_klf, (7.20) 



we have defined K — ^K±K\\ and from Eq. 5.63| , the 
three-dimensional logarithmic coupling constant is 



ao 



1-fao ln(A/fc_L)' 



(7.21) 



Eqs. 



7.1£ 



7.21 



7.1C 



The 
time- 



have simple physical significance. The 
first term in each of the square brackets in Eq, 
represents time-dependent "thermal" fluctuations, 
second, (5-function term represents a static, 
independent distortion of the smectic, and is in fact iden- 
tical in form to the perturbative expressions of Section 
IV. This is roughly because at the uv momentum cut- 
off (where the rescaled correlator is evaluated), the dis- 
placement fluctuations are strongly suppressed, and dom- 
inated by the value of the random force at origin = 
(as in the naive Larkin approximation). 



CEA(r) 



A 



(d-2)/4 / I 1^1/2 



1/2 



(7.23) 



where the scaling function obeys {x 0) constant 

and f'^\x oo) — > x^^'''. For d < 3, in the asymptotic 
limit x > K^^/^v, C^{r) becomes infinitely anisotropic 
owing to the difference between the convective behavior 
along the direction of motion (x) and diffusive transport 
transverse to it (along rj_). We find 



CEA(r) = A 



3-d 



"fvK 



„(2 



A 



(7.24) 



where the asymptotics of the scaling function are 



x) (X 



constant, for x 







(7.25) 



In three dimensions, the power-law in Eq. 7.24 should be 
replaced by the logarithm {r±y^ In |Arj_|. 

The other important physical correlator measures the 
thermal fluctuations around the static distortion mea- 
sured by Cea- The thermal correlator is naturally de- 
fined as 



CT(r,i) = C(r,i)-CEA(r). 



(7.26) 



For simplicity, we consider only the equal-time thermal 
fiuctuations; corresponding results for non-equal times 
are easily obtained from Eqs. 7.1E-7.21. In two dimen- 



sions, the thermal correlator is logarithmic. 



CT(r,O)~2reffln0b 



(7.27) 



where K = ^Kj_/K\\. Note that the (in 2d nonuniver- 
sal) coefficient of the logarithm above is proportional to 
the effective temperature, and therefore bounded below! 
This is a consequence of the roughly semicircular RG 
fiows in Fig. 6. It is interesting to further note that by 
extending above calculations to 2 + e dimensions, using 
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results of Sec.VI, we find that the logarithmic growth 
of thermal correlation function is infact super-universal 
(i.e., independent of d) for 2 < d < 3, with the prefactor 
Tgg above replaced by a universal fixed point value of 

T* = 1 - 2e (Eqjel^). 

In thr ee dimensions, the logarithmic decay of a in 
Eq. 7.21 re nders the Fourier transform of the thermal 
term in Eq. 7.19 non-divergent. The thermal correlator 



thus saturates at long distances in d = 3, 



CT(r, 0) |r| — oo Ct,0 oc 



ao 



Inaol 



(7.28) 



Comparing the static (EA) and thermal correlators 
above, we see that at any fixed time, the static contribu- 
tion to C(r, t) dominates at long distances, i.e. 



C(r,t) irHoo CEA(r)- 



(7.29) 



Thus the naive perturbative results of section IV are es- 
sentially correct for the long-distance properties of the 
correlation functions. In particular, the structure func- 
tion in the smectic phase displays power-law smectic 
Bragg peaks (translational QLRO ) in three dimensions, 
with fully rounded smectic peaks (translation SRO) in 
two dimensions. 

Similar analysis of the response function -R(k, uj) shows 
that to the leading order in disorder it is given by 



i?(k,cj) = 



1 



(7.30) 



and is identical to the response function in the linearized 
theory of Sec. HI, but with a drag coefficient 7 multiply- 



ing u) (but not kx) replaced by the disorder enhanced 
renormalized 7_R(fcj_), 



lR{k±) = 76 



*(A/fc_L) 



(7.31) 
(7.32) 



Since 77?, is finite, the renormalized response function 
i?(k, Lo) implies an analytic response to a uniform trans- 
verse external force, and a finite linear mobility 7^^ in 
the limit of vanishing transverse force. 

We now turn our attention to the nonlinear dynamic 
response to a transverse force f±. A true calculation of 
the non-linear response to f± is beyond the current ca- 
pabilities of these methods. A relatively simple scaling 
argument suffices, however, to obtain a rough picture of 
the dynamics. Consider introducing an explicit force into 
the MSR functional. Since it couples directly to 0, sim- 
ple rescaling and use of the relations between the scaling 
exponents X: ^ C leads to a strong growth of the 
rescaled force f±{b) under the RG, 



(7.33) 



When this rescaled force becomes large, of the order of 
the typical depinning force fc ^ A''/^-\/A"(0), it be- 
comes a strong perturbation and the damping coefficient 



7 should cease to renormalize. Choosing = fc de- 
fines a (/j_-dependent) rescaling factor b at which to eval- 
uate 7 via Eqs. 7.10- [7.li . This gives an effective drag 
coefficient 

7cff(/±) = 7fle-''(/^//=) , for < /, , (7.34) 
= 7, for/i>/c, (7.35) 

where v = aK^^K^A"^ / {•jv}'^ . Combining this result with 
the definition of 7 



v±if±) 



fi 



7off(/-i 



(7.36) 



we obtain a nonlinear response which interpolates be- 
tween two different linear responses at small and large 
/_!, schematically displayed in Fig. 7. We remark in pass- 
ing that this simple argument neglects possible renormal- 
izations of f± by the random force. This occurs at zero 
temperature, and allows for the possibility of a true crit- 
ical transverse force in that case. A careful treatment of 
the finite-temperature crossover to the above (generic) 
behavior is still an open problem. 




Fig. 7: Schematic of the nonlinear transverse re- 
sponse for the moving smectic (neglecting the permeation 
mode). Including the permeation mode simply adds a 
linear function of slope Uperm. to the above plot. 

The finite-temperature "IV" (velocity-force) curve 
gives an interpolation between two linear regimes, above 
and below fc. For small 7W, the renormalized mobility 
7^^ is quite small. In this case the IV has a non-linear 
feature similar to that of a threshold, but nevertheless 
completely analytic. In addition, the permeation mode 
provides a second channel for transport, which further 
enhances the linear mobility. A possible (though highly 
speculative) scenario at zero temperature is a small linear 
mobility via the permeation of a small concentration of 
fluctuating defects "activated" through local chaotic dy- 
namics, superimposed upon a sharp threshold for motion 
of the smectic density wave. However, we again caution 
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(as in the Introduction), that zero temperature dynamics 
may be highly nonuniversal, with quahtatively different 
behavior occurring in different systems. 



VIII. DISCUSSION 

A. Relation to previous work 

Much of the recent interest in driven solids stems from 
the work of Koshelev and Vinokur (KV)J23 who applied 
perturbation theory to compute the mean-square dis- 
placement of the driven solid and argued that the driven 
system might exhibit a nonequilibrium phase trapsition 
from a moving liquid to an ordered moving solid. Ej The 
argument was based on the notion that at large velocities 
the effect of pinning could be described by an effective 
"shaking" temperature Tgh ^ With this assump- 

tion, KV recast the problem into equilibrium form. This 
suggested only thermal roughness for the driven lattice, 
so that the mean square displacement would always be 
bounded in c? > 2 (with translational QLRO in d = 2), 
arguing indeed for the stability of the driven lattice for 
d>2. 

Experiments on flux lattices in type-II superconductors 
have indeed shown evidence for current-induced ordering 
of the vortex array. This evidence has beeH-obtained both 
indirectly through transport experimentsEJ-ell and by di- 
rectly probing the stajcture of a driven vortex lattice 
by neutron scatteringuS and decoration experimentsc3'E3. 
Numerical simulations of driven vortex arrays in two di- 
mensions also provide clear ipdicatJpa of ordering of the 
sliding lattice at large drivesCJO'Eao. We discuss these 
experiments in more detail below. 

The notion that non-equilibrium effects might plajtan 
important role was addressed by Balents and Fishenlj in 
the simpler context of CDWs. They classified the pos- 
sible phases of driven CDWs and showed that the KV 
predictions were violated in this case. In fact, certain 
static components of the quenched disorder persist in a 
coarse-grained model even at very large velocity. This 
static random force arises physically from spatial inho- 
mogeneities in the impurity distribution and represents 
a sort of random drag. Therefore noise in the nonequi- 
librium steady state of the sliding CDW never mimics 
thermal noise, which is uncorrelated in time. As a re- 
sult these authors concluded that the analogous moving 
solid phase of a CDW is stable in 3d at large velocities, 
but does not exhibit the true long-ranged translational 
order of an equilibrium 3d crystal. Rather it exhibits the 
algebraic decay of correlations characteristic of a 2d equi- 
librium crystal (QLRO). In rf = 2 in contrast the moving 
CDW appears unstable to the proliferation of phase slips. 

A similar reinvestigation of the driven lattice was ca- 
cently undertaken by Giamarchi and Le Doussal (GL)E3. 
These authors focused on the behavior of fluctuations in 
the direction transverse to that of the driving force. They 



pointed out that periodic static component of the pinning 
force persist in the transverse direction at large drives 
and suggested that motion in the driven solid occurs 
along elastically coupled channels parallel to the direction 
of motion. These authors assumed that disorder-induced 
displacements in the direction parallel to the driving force 
always remain bounded and neglected all fluctuations in 
this direction. They also suggested that the signature 
of this anisotropic sliding state, the anisotropic moving 
Bragg Glass (BG), would be the existence of a finite 
transverse critical force at zero temperature. 

Considerable elements of GL's work are retained in our 
calculations. The static periodic transverse force is also 
non- vanishing in our model, and leads to glassy channel- 
like motion at very low temperatures. There are, how- 
ever, several key differences from GLs work. First, we 
give a model independent characterization of the possi- 
ble phases. Second, we derive our equations of motion 
from a microscopic dynamics. These have a truly non- 
equilibrium form, which cannot be obtained simply by 
a Galilean boost of the equilibrium equations of motion. 
Third, we argue that the most stable driven phase is the 
transverse smectic, with short-range rather than long- 
range longitudinal order. Fourth, we point out the im- 
portance of the permeation mode in the transverse smec- 
tic, which implies a non-zero transverse linear response 
at any finite temperature. Lastly, we show that the trans- 
verse displacements themselves are much less glassy than 
in equilibrium, owing to the strong enhancement of ther- 
mal noise caused by the breakdown of the fluctuation 
dissipation theorem. 

Recent simulations in d = 2 have confirmedi— thei 
anisotropic channel structure of the sliding stateBa'cSl 
These simulations also indicate that dislocations with 
Burger's vectors parallel to the flow are unbound, so that 
the intra-channel order is liquid-like. The driven array 
thus indeed appears more consistent with a transverse 
smectic phase than the moving Bragg glass. 

All the aforementioned analytical treatments n eglec t 
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KPZ-type nonlinearities (C coefficients in Eq, 
These are perturbatively irrelevant in the smectic state, 
as argued above. They probably play a role in the be- 
havior of longitudinal correlations, and may well deter- 
mine the length scale at which longitudinal dislocations 
appear in the moving lattice, i.e., the scale at which a 
driven lattice is unstable to the transverse smectic. Such 
effects have recently been studied in a simpler model rel- 
evant for CDW motion, which contains only longitudinal 
degrees of freedom.c3 



B. Experiments 

As discussed in the Introduction, there are many phys- 
ical realizations of dirty driven periodic media. Among 
these, the magnetic flux lattice in type-II superconduc- 
tors is perhaps the system that has been most studied 
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experimentally in recent years and where our predictions 
can most easily be tested. 

The large majority of experimental work has focused 
on the nonlinear transport properties of these systems 
for driving forces near the zero-temperature depinning 
threshold. Our work, in contrast, focuses on the proper- 
ties of the sliding state well above threshold, where the 
velocity-force characteristic approaches a linear form. In 
this regime one is rather interested in the positional and 
temporal order of the sliding medium. 

The positional correlations in a current-driven mag- 
netic flux-lattice can be studied directly both by numeri- 
cal simulations and experiments. Numerical work is very 
useful as it can provide both direct real space images of 
the sliding lattice as well as quantitative structural infor- 
mation like the structure factor. Recent simulations of 
two-diinensional flux lattices {dt = 2, di — 0) by Moon 
et al.ES are in agreement with our finding that the pe- 
riodicity of the driven flux lattice along the direction of 
motion is absent, i.e., the sliding lattice is a smectic. Real 
space images of the driven lattice show that motion oc- 
curs along channels that are aligned with the direction 
of the driving force and periodically spaced in the trans- 
verse direction. Phase slips, however, occur at the chan- 
nel boundaries, indicating that the channels are uncor- 
related and the longitudinal structure is liquid-like. The 
structure factors obtained from these simulations show 
sharp algebraically-divergent peaks at the reciprocal lat- 
tice vectors Q normal to the external drive. The peaks at 
the other reciprocal lattice vectors have a very small in- 
tensity that decays exponentially with system size. ^'Y)'". 
lar results have also been obtained by other authors 
One detailed aspect of the results of Ref. ^ which docs 
not agree with our theoretical expectations is the alge- 
braic decay of the smectic Bragg peaks. Our theory pre- 
dicts in fact (possibly stretched) expon entia l decay, due 
to the linear displacement growth in Eq. 7.24 . We believe 
the observed power-law structure factor scaling is likely 
a crossover phenomenon, perhaps enhanced by the dis- 
persive elastic moduli due to the short-scale logarithmic 
character of the inter-vortex interaction. 

One of the first direct experimental evidence of the or- 
dering of the sliding fiux lattice at large velocities was ob- 
tained some time ago by neutron diffractionE^ These ex- 
periments have not, however, been able to quantitatively 
determine the structural properties of the driven state. 
More recently, the channel structure of the driven flux 
lattice was observed directly by decorating the current- 
driven fiux array in iV6>5'e2.c3 By digitizing the decoration 
images Pardo et al. have also very recently obtained the 
structure factor of the driven array that indeed has the 
transverse peaks characteristic of a smecticEj. 

Temporal order in the driven medium should manifest 
itself in narrow band noise (NBN) and mode-locking phe- 
nomena. Both of these have been studied extensively in 
the context of charge density waves, but have not been 
observed in flux lattices, indicating that the driven flux 
lattice lacks long-range temporal order. The spectrum 



of voltage fluctuations can be probed by applying a dc 
current / that yields a driving force / ^ J on the flux 
lattice. The local field induced by flux motion is given 
by 



Q 



(8.1) 



with (pQ the fiux quantum. In a moving solid we expect 
(p — vt + 4>. In a perfect lattice the local voltage con- 
tains therefore oscillatory components at the frequencies 
luq — Q-v. The Fourier spectrum of voltage fluctuations, 
defined as, 



Sico) 



(8.2) 



will then contain, in addition to a dc component, a sharp 
fundamental peak at ui — Qi • v, with Qi ^ 27r/a, and 
smaller peaks at all the higher harmonics. This type of 
spectrum, usually referred to as NBN, has indeed been 
observed in sliding CDWs and is generally regarded as 
the signature fOf the absence of appreciable phase slips 
in the system.EJ In contrast, in a liquid phase we expect 
the Fourier spectrum to have a Lorentzian form, centered 
at the frequency uji. A broad power spectrum of volt- 
age fluctuations, with an uj~°' decay at large frequencies, 
known as broad band noise (BBN), is observed in CDWs 
when macroscopic velocity inhomogeneities p:ising from 
phase slips are present in the driven systemO NBN has 
not been observed in current-driven fiux lattices. This 
is consistent with our finding that correlations along the 
direction of motion are always liquid-like, indicating lon- 
gitudinal phase slips are present in the system. 

Systems with temporal LRO should also exhibit com- 
plete mode-locking. Mode locking is an interference effect 
that can occur in a periodic medium driven by both a dc 
and an ac force. Keeping fixed the amplitude and the fre- 
quency of the ac force and varying the amplitude of the dc 
component, one observes steps in the dc response, known 
as Shapiro steps. These steps arise from mode locking of 
the frequency of the applied ac force with the internal 
oscillation frequencies loq of the periodic medium. Com- 
plete mode locking (steps in the dc response that remain 
constant over some finite range of dc bias) to an arbi- 
trarily weak external ac drive has again been observed 
in CDWs.E3 Assuming, as our calculations suggest, that 
the driven flux lattice has only short-range longitudinal 
order (i.e. SRO for wavevectors with non-zero wq), we 
would expect at best incomplete mode locking above a 
non-zero (perhaps large) threshold ac drive. For (2 H- 1)- 
dimensional flux lattices again no complete mode locking 
is expected in the smectic state, as the driven lattice only 
has quasi-Bragg peaks along the direction in wavevector 
space perpendicular to the velocity. If a longitudinally 
ordered phase ( "Bragg glass" ) were stable in 3d, it ought 
to exhibit mode locking and NBN; the available experi- 
mental evidence seems not to support this possibility. 
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Finally, it has been suggested that a sliding flux lat- 
tice will exhibit a finite threshold force for response to 
an addiiional driving force f± transverse to the mean 
velocityE£lE3 and no transverse linear response at zero 
temperature. The behavior in a purely dissipative (over- 
damped) model at zero temperature is somewhat non- 
universal, and indeed, such a transverse critical force is 
certainly likely in many possible phases, including both 
the smectic and Bragg glass (BG). At finite temperature, 
a sharper distinction can be drawn. A naive extrapola- 
tion of the Bragg glass theory to finite temperature would 
predict an exponentially small, non-linear transverse re- 
sponse at finite temperatures. 



,,BG 



exp 



(8.3) 



with/n" cx 1/T at low temperatures and small f± 



Eq. 8.3 relies, however, upon two assumptions: (1) den- 
sity fluctuations (the permeation mode) can be neglected 
and (2) the T = fixed point is stable. Our works shows 
that both these assumptions are incorrect and invalidates 
Eq. ^.3| . As discussed, in Section VI, we then expect a 
velocity-force curve with the crossover behavior shown 
in Fig. 7. For low transverse driving forces there are two 
small (but nonzero) linear components to the mobility, 
i.e. 



(McIw. + Mpcrin.(^)) f±- 



(8.4) 



The mobility of the density wave, fidw., has been esti- 
mated in Sec. VII. It is activated at low temperatures, 
and also suppressed for small sliding (longitudinal) ve- 
locities, with the form 



In' 



^-lg-<TK|,if_LAV(71')' 



(8.5) 

where a cx 1/T at very low temperatur es (m ore details 
are given in the discussion following Eq. 7.12). The sec- 
ond term is a mobility associated with the permeation 
mode: 



A'pcrm. 

(T) ^ floe 



(8.6) 



for T <C Ud- Although we have not analytically de- 
rived the above exponential form of the mobility, it 
seems extremely likely on physical grounds. Such be- 
havior arises from an activated concentration of mobile 
defects (vacancies or interstitials in 2d, kinks and/or va- 
cancy/interstitial lines in 3d) with finite "energy" cost 
Ud- These flow linearly in response to a driving force with 
mobility fiQ. Additional non- linear behavior will be su- 
perimposed upon this linear term, but it is subdominant 
at small f±. At low temperatures these non-linearities 
sharpen to a threshold-like feature around a finite f±c, 
so that the identification of an experimental system as 
either a BG or a smectic by the transverse "threshold" 
or critical current may be misleading. 



C. Open questions 



We conclude with a summar y of s ome of the many re- 
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provide for the first 



maining open questions. Eqs. 
time a complete set of hydrodynamic equations to de- 
scribe nonequilibrium states of driven periodic media. A 
systematic RG (or other) treatment of them is, however, 
daunting. Aside from the obvious algebraic complexity, 
the structure of the equations also raises some interesting 
conceptual issues. Naively, since first derivative terms are 
present in all the variables, all spatial coordinates seem 
to scale isotropically, with dynamical exponent z = 1. 
This power counting, however, is inconsistent with di- 
rect perturbative (and more sophisticated) calculations 
of physical quantities. Preliminary investigations of this 
problem suggests novel scaling without unique anisotropy 
or dynamical exponentsEj A more modest goal, requir- 
ing only straightforward (if rather tedious) calculations 
is to extend the RG treatment of the smectic to include 
the permeation mode, i.e., th e couplin g to density fluc- 
tuations, as described by Eqs. 4.6-4.7. 

A systematic treatment of topological defects in dis- 
ordered lattices is also lacking. Some recent progress 
along these lines was made recently for equilibrium elas- 
tic glasses in Ref. ^ - it remains to be seen whether this 
work can be extended to driven systems. When disloca- 
tions are present j_inany interesting questions remain at 
zero temperatureEj. A possible plastic depinning tran- 
sition and its properties - the associated scaling, noise, 
and hysteresis - are not presently understood. 

Finally, an extremely interesting line of inquiry is to 
explore how some of the techniques used here might be 
extended to inertial models appropriate for friction and 
lubrication. Such systems presumably are controlled not 
only by the physics described here, but also by nontrivial 
mechanisms of dissipation and possibly chaotic dynam- 
ics. 
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APPENDIX A 



In this appendix we carry out the mode-elimination 
needed to coarse-grain the driven lattice in a periodic 
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pot ential . Expanding the e '^^ within the brackets in 
Eq. 2.22| to first order gives 



u<i{r, t) +u>i(r,t)] 
X i^, [x + u< (r, t) + u> (r , <) + vt , z] 



(8.7) 



o> 



This average can be evaluated by expanding the force F 
in u> and averaging over the fast modes term by term. 
However, in the absence of thermal noise, the only non- 
zero fast correlator is the response function, which van- 
ishes by causality at equal times. Since all the fields in 
the first order term are at equal time, all the terms in- 
volving any u> or u> fields are zero, and this average 
simply gives 



half of the square of 5S),J is su btracted off, as demanded 
by the logarithm in Eq. ( 2.22| ). As we saw earlier, only 
non-equal-time response functions can survive the aver- 
age. This can occur only via the first order expansion of 
the displacement field out of one of the force terms, to 
be contracted against the response field at the non-equal 
time. Taking into account the two ways in which this can 
be achieved gives the result 



^c.(2) _ 
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xx' 



u< 



- vt, z] 



xF_,-[x' + u'< +vi',z']Gjfc(r-r',i-t'). 
Here the response function is (in matrix form) 
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off. 



^( {t<,(r,i)#,[x-hu<(r,<) +vt,z]. (8.8) G(q,t^) = 



n -1 



iujl + K{q) 



e{\q\-K). 
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To leading order, then, the force is unrcnormalized. 

The first non-trivial correction arises at second order. 
This effect is actually physically transparent. At first 
order the fast modes simply respond linearly to the adi- 
abatic motion of the slow modes. However, because of 
the mode coupling, this response is fed back and felt at 
second order by the slow modes again. It is this feedback 
that corrects the motion of the slow degrees of freedom, 
and is perturbatively estimated in the next correction. 
This is given by 



SS. 



(2) 
eff. 



The ^-function indicates that only a partial mode elim- 
ination has been performed, so that the slow modes re- 
main as dynamical variables in the coarse-grained the- 
ory. Generally, however, the expressions obtained in this 
section have smooth limits as A — > 0, and are well ap- 
proximated therefore by using the response function in 
the full Brillouin zone. 

Using the Fourier decomposition, Eq. 2.23, the second 
order contribution to the effective action can then be cast 
into the form 



SS. 



(2) 



E 



xFi[x + u + vt, z]Fj [x' + u' + vt', z'] 
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(8.9) 
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where we have indicated the arguments of u and u by 
the presence or absence of a prime, and furthermore sup- 
pressed the mode decomposition. The superscript c indi- 
cates a cumulant, or disconnected average, meaning that 



where Sfi may be interpreted as an additional effective 
force in the coarse-grained equation of motion for the 
slow modes u<. The dominant terms in Sfi are those 
which do not oscillate as x is varied. Keeping only these, 
it takes the form 



<5/.[u<,r,<] = ^ f ^zQ,Q,Q,e'Q(--'+"-"'+-(t-*'))|C/Q|2G,.fc(r-r',i-0. 

X' Q 

This is greatly simplified by gradient expanding the difference in displacement fields 

u - u' w (r - rTdaU + {t - t')dtu - i(r - r')"(r - r'fdc^dpu. 
We then obtain an expansion for the force corrections, 

= E / E*^?'Q^^?'«if^Qi'^''^'^'''^"*'^G^fe(^''^')-^[u<], 

X' Q 
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(8.15) 



where 



^■[u] 1 + iQi 



-\QiQ,nr'''rPdc.u'dpu^ + ■■■. 



(8.16) 



Including these corrections into the effective action gives Eq. 2.24 in the main text. 
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APPENDIX B 

Here we present the details of niode elimination for 
a lattice driven over a disordered potential. Expanding 
from the exponential (e~'^) and averaging it pertur- 
batively over the modes outside of the cut-off gives, to 
hncar order in 5*1, 



(5: 



6S. 



(la) 



ss. 



(lb) 

:ff. ' 



where the first term is, as before 
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and simply returns the uncorrected bare random force. 
The next correction is 



SS, 



(16) 
eff. 
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u\r,t)Gjk{r-r',t-t') 



xdidjdkt[r - r' + u(r,t) - u(r',t') + v(t - t')]. (8.19) 



Using Eq. (8.14), this becomes 



x,x' 



fyfei [x' - voi',z'] 



(8.20) 



r "daUi + t'dtui - \r ^dadpui 



where we have abbreviated F 



didj ■ • • r. This cor- 



rection to the effective action again represents gradient 
terms in the coar se-grained equation of motion, of the 
form of Eq. ( ^.24| ). Extracting these coefficients, we ob- 
tain the formulae quoted in Section 11. B. 2 for the deriva- 
tive coefficients 7, A, B and C. 

It remains to consider the renormalization of the force 
term itself. This vanished in the case of the periodic 
force, due to the (assumed) incommensurability of the 
lattice and the substrate. The random potential, how- 
ever, has Fourier components commensurate with the 
moving lattice, which thereby generates such a renormal- 
ization. To evaluate it, however, we must go to higher 
order in the disorder variance F. In particular, we con- 
sider the second correction term 



SS, 



(2) 



.21) 



This is explicitly 




SSf^ = y ] { MiU^Ftj[xi-X2 + Ui-U2+v(ti-^2),Zl-Z2] 



X'U3'U4Ffc;[x3-X4 + U3-U4 + v(^3 - t4),Z3-Z4] 



.22) 



0> 



where we introduced the obvious abbreviation for the four lattice sums and longitudinal space and time integrals. 
This expectation value contains several terms, depending upon the number of u fields which are contracted to give 
response functions. Terms with one contraction leave 3 response fields, which represents the generation of a skewness 
to the distribution of the random force, and can be neglected in what follows. Forming three contractions leaves a 
single response field, which will give higher-order corrections to the coefficients determined above, and can thus also 
be neglected (for weak disorder). Forming two contractions leaves two response fields, which is of the proper form to 
renormalize the random force. 

To obtain these terms, the F's must be expanded to second order in the fast fields u>. This gives 



SS. 



(2) 
off. 



(12)f H„(34) (u\uiiuTy - «">)ii34«> - u 



+F,,™„(12)Ffc,(34) (ul4«> - )«> - u^y)u^uiy 



4>, 



>0 



/ >0 



(8.23) 



where we have now also abbreviated the arguments of the force correlators. The contractions inside the angular 
brackets can still be performed in several ways. Each such choice gives rise to a separate term containing two response 
functions and a combination of slow fields at different space-time points. To determine the desired correction to the 
random force correlator, we keep only the leading term in a gradient expansion of the slow fields (i.e. zeroth order in 
the gradients). Lengthy but simple calculation gives 



^^cff. = / ^iXr,t)^^(r',0'5ry[x-x' + u<-<+v(i-i'),z-z'], 

^ ^ ^, Jzz'tt' 



(8.24) 



where the renormalization of the force-force correlator is 



31 



Sfij [x, z] 



>^r,:fe„i[x-Xi-vti,z-zi]rj7„[x+X2+vi2,z+ Z2]G'fe„(l)G'/„i(2) 

"12 

+ 2Vfifcm[x-Xi-Vii,Z--Zi]fj7„[x2+vi2,Z2]G/cn(l)G/m(2 - 1) 
"12 



+f jj>„„[x,z] y^ffei[xi-X2+v(ti-t2),Zl-Z2]Gfcm(l)G;„(2) 



-Xi-X2 + v(ti-t2),Z + Zl-Z2]Gfe™(l)G,„(2) 



f [Xi+ Vti , Zijffej [x+X2+ Vt2 , Z+ Z2] G;,„ (1) Gfc„ (2) 



•'12 

- f i;„i„ [xi+ vi 1 , zi ] f [x+ X2+ vi 2 , z+ Z2 ] G/,„ ( 1 ) Gfc„ ( 2 - 1 ) . 
Integrating this horrendous equation over r, one obtains the formula for g^j quoted in the text. 



(8.25) 



APPENDIX C 

Here we estimate the correlator of the nonequilibrium 



part of the static random force, gij, given in Eq. 2.51 



ffi = / (q* • v)^gfeg/9„i9„|r(q)pGfo„(q, qt • v) 

X [G,„(q, -qt • v) - G,„(q, q* • v)] . (8.28) 



the limit of large sliding velocity v. First we note that as 
gij is a symmetric tensor, it can be written as 



where v = v/w and 



.26) 



50 + .91 = / q gfeg/gmgn|r(q)| Gfo„(q, q* -v) 



[G/„(q, -qt • v) - G;„(q, qt • v)] , 



and 



The physical case of interest here is the one where the 
range ^ of the pinning potential is small compared to the 
lattice constant, a. The Fourier transform of the vari- 
ance of the random potential, r(q) is then a very broad 
function on the scale of the Brillouin zone. To proceed, 
we exploit the periodicity of the elastic propagators in 
reciprocal space, by letting c\t = k + Q, where k only 
spans the first Brillouin zone. The integral over qt is 
then replaced by an integral over the first Brillouin zone 
27) and a sum over all reciprocal lattice vectors. Focusing for 
concreteness on the evaluation of gi {go can be evaluated 
by a similar procedure), we obtain 



9i = Y. I I [(k + Q)-v]'(fc + Q)fc(fc + Q);(fc + Q)„(fc + Q)„|r(k + Q)|2 

cGfc™(k, g,, (k + Q) • v) [G,„(k, g,, -(k + Q) • v) - G,„(k, g„ (k + Q) • v) 



(8.29) 



We now split the reciprocal lattice vector sum in Eq. ( g.29 ) in two parts by separating out the terms with Q • v = 
and write 



.91 - .9i + 9i ■ 



(8.30) 



The term g\ denotes the contribution to gi from the sum over reciprocal lattice vectors satisfying Q • v 7^ 0. In this 
term we neglect everywhere k compared to Q and obtain 

9i'^^ J2 (Q-v)'QfcQiQmQ„|r(Q)p /" Gfc„(k,g,,Q • v) 

Qv#0 "'''■«- 



[G,„(k,9.,-Q-v)-Gi„(k,g„Q.v)]. 



5.31) 
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In the limit v >> 2TT^/{ac), with c a typical elastic 
constant, we ca n now approximate the elastic propaga- 
tors in Eq. ( ^.31 ) by neglecting the in-plane elastic matrix 
compared to the frequency 7Q • v, i.e., 



Gfo„(k, g^, Q • v) 



C44g; 



i-fQ ■ 



..32) 



By inserting Eq. (p.32|) into Eq. (8.31), we obtain 



9i 



(1) 



2po 



E 

Q-v#0 



(Q-vo)V|r(Q)|^ 



(7Q 



[M.r + (7Q • v)2] 



212 ■ 



.33) 



For di = this gives 



2po 

ry2y2 



Q-V05^0 



0'|r(Q)|^ 



/ r(Q = o) y 1 

V 7u / 



(-)' 



(8.34) 



where A w r(Q = 0)a/^ is the variance of the equilib- 
rium part of the static pinning force defined in Eq. 3.8 . 
For di = 1 the gz-integral is easily performed, with the 
result, 



3/00 



1 



IQxlQlr(Q)|= 



A2 



(a7t>)3/2^C44a$2 



(2) 

In the contribution g]^ ^ from the Q • v = of the 
reciprocal lattice vector sum the integral over k is domi- 
nated by k near the center of the Brillouin zone. In this 
term we can therefore approximate the elastic propaga- 
tors by their long wavelength form, given in Eqs.|3.13|and 
3.14 . After some lengthy algebra, one can show that in 
the limit of large shding velocities, v » 2TTj/{acQQ), 

~ (^/<^)ffi^^ is therefore negligible for short- 
ranged pinning potentials. 

The evaluation of go can be performed by the same 
method, with the result 



.35) 



90 
90 



91 



d/ =0, 



5iln(a/0, di = l. 



In summary, we find that at large sliding velocities 
both components go and gi of the correlator of the 
nonequilibrium part of the static random force have 
the same asymptotic dependence on v and the disorder 
strength A, with 



9o,i 



A2 



(2dt-d,)/2- 



V 



..37) 



APPENDIX D 



Here we examine the predictions of the perturbation 
theory described in Section III for the real-space decay 
of positional correlations in a dt-dimensional lattice of 
magnetic flux lines {di = 1), driven in the x direction. 
To obtain the mean-square displacement in real space, 
we need to evaluate integrals of the form. 



B{r) - 2 



1 — cos(q • r) 



(8.38) 



where c stands for either cge or cn -I- cge and the prime 
denotes a cutoff at = A. The integrals over and 
qz are easily performed. Letting u = q±x_\_, one obtains. 



B{v) 



\y 



\2-dt 



.39) 



where s = \x\c / {vjy'^) and C ~ ^ c/ C/^/i\z\/\y\. The scal- 
ing function J-['^*\s,C,,K\y\) is given by 



.f(*)(s,C,A|x^| 



dvL 



,dt-3 



osh(c.)-Vc«<i>(Vi.-^) 



2Vi 



1 — cos(y • u) 



(8.40) 



Here u is a {dt — l)-dimensional vector, with u = u/m, 
and $(a;) is the error function. We are interested in the 
asymptotic behavior of the scaling function for dt = 2. 
For C = we find 



J-f)(s-.0,0,A|2/|)~ln(A|y|), 



^f)(.»l,0,AM)^ln(^:^) 



.41) 



(8.36) For C >> 1, or |z| >> 12/1^044/0, 



^P'(s»l,0,A|,|)~ln(ii|). 



..42) 



The scaling of the mean square displacement is therefore 
anisotropic, but logarithmic in all directions. 
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APPENDIX E 

In this appendix we outline the details of the 2 + e 
RG calculation for the single Fourier mode driven smec- 
tic model defined by Eq.6.7 in Sec. VI of the main text. 



It is convenient to employ the Martin-Siggia-Rose (MSR) 
formalisms. In this formalism one studies the dynamic 
generating functional Z which is a trace over the displace- 
ments (^(r,i), with the constraint that 0(r, i) satisfies 
the equation of motion Eq j6.7| , imposed via a functional 
(5- function as an integral over a response field (/)(r, t). Av- 
eraging over the noise r/(r, t) and the quenched random 
force Fp[0(r, i), r], the problem can be recast in the form 
of a dynamical field theory, 



-So[0,0]-Si [0,0] 



(8.43) 



where in addition to the standard quadratic part of the 
action 5*0 



Sq — 



(r, t) {jidt + vd,) - {K\\dl + K^^D] 0(r, t) 



-lT4^{v,tf , (8.44) 
there is a contribution Si due to disorder 
Si^-\ I 0(r, t)^{r, t')Ai cos[go(0(r, t) - 0(r, t'))] . 

^ Jr,t.,t' 

(8.45) 

In the above, after averaging over disorder we have kept 
only the most relevant lowest Fourier component of the 
random force correlation function. We have also de- 
formed the functional integral contour over the response 
field (j) to the imaginary axis. 

We employ the standardjnomentum shell renormaliza- 
tion group transformationEJ, by writing the displacement 
field as 0(r, t) — cj)^ (r, t) + (p^ (r, t), integrating perturba- 
tively in Ai the high wavevector field (/>'*(r, t) nonvanish- 
ing for Ae~' < gj_ < A (with no cutoff on the momentum 
along the direction of motion and on w), and rescaling 
the lengths, the time, and long wavelength part of the 
fields with 



i)<{v,t) = eX' 



0<(r,i) = e 



(8.46a) 
(8.46b) 
(8.46c) 
(8.46d) 
(8.46e) 



so as to restore the ultraviolet cutoff back to A. Because 
the random-force correlator Ai term is a periodic func- 
tion of (j), it is convenient (but not necessary) to take the 
arbitrary field dimension x = 0, thereby preserving the 



period 2Tr/qo under the renormalization group transfor- 
mation. Under this transformation the resulting effective 
free energ y function al can be restored into its original 
form Eqs.B.44-B.45 with the effective Z-dependent cou- 



plings. To zeroth order we obtain: 
d7(0 



dKj_{l) 

di 



-(rf± + C + x)7(0 

= {dj_ + z + xhv{l) 
{d^+C + z-2 + x)K±{l) 



dl 

d'yv{l) 
dl 



dK\\{l) 
dl 

dAijl) _ 
dl ' 
dTjjl) 
dl 



= (d± - c 

{d±+C + 2z 



+ Z + 



X)K\\{1) 
2x)Ai(0 



= {d^+C + z + 2x)T^{l) 



;.47a) 
.47b) 
!.47c) 
.47d) 
;.47e) 
5.47f) 



As in the calculation of Sec.^ statistical symmetry un- 
der an arbitrary time-independent shift of the displace- 
ment field 0(r, t) ip{r, t) + /(r) (time-translational in- 
variance) guarantees that 7D, K± and i^y do not acquire 
any graphical corrections, i.e. their flow equations above 
are exact. Imposing this requirement at the tr ee level on 
the first two coefficients, and using Eq. 8.47a we obtain 
C = 2 and x = —d± — z. 

More generally but equivalently we look at the dimen- 
sionless coupling constants 



T 



2(Xiif||)i/ 



d-2 



Ai 



d-3 



K±"f\v\ 

which have tree-level fiow equations 
f = (2-d)T(0 



dAi 
dl 



(3-d)Ai(0 , 



(8.48) 
(8.49) 

(8.50) 
(8.51) 



whose flow is independent of the arbitrary choice of 
rescaling exponents appearing in Eqs. j.46a- ^.464 

We now proceed to higher order in Ai, perturbatively 
integrating the short length modes (j)^{r,t) and (j}-^{r,t) 



[d<l)<d(t> 



= / [d^<d(b<] e-5o[0<.0<]-5S[0<,0<] ^ 



(8.52) 
(8.53) 



where the graphical correction SS to the action (dropping 
an unimportant constant) is 



,55[(^<,0<] = (5l[0,0]); 



2('5'i[0,' 



(8.54) 
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where the subscript "c" means cumulant average, and 
the averages are performed with the quadratic action 
with correlation and response functions, C{q,uj)V = 
((/>(q,w)(/)(-q,-w)), G{(i,uj)V = ((?!)(q,tj)0(-q, -w)), 
respectively given by 



C{q,u;) 



2T7 



(8.55) 



which can be read off fro m E q j8.44 . Although naively 
one would expect from Eq. |6.6| that as in the equilibrium 
problem for d > 2 temperature is an irrelevant variable, 
as we demonstrate below (con sistent with the functional 
RG treatment of the Sec. VP this zero temperature fixed 
point is destabilized by the finite velocity motion. We 
will therefore work at a finite temperature. 

The first order correction {Si[(j), 4>])y contributes to the 
renormalization of 7, T7 and Ai, which we designate as 
5^^'^\ S(Tj)^^\ and (JA^^-* and illustrate graphically in 
Fig. 





Fig. 8: Two diagrams that contribute to the renormal- 
ization of Ai, T7, and 7. The full line corresponds to the 
correlator C^, the full- wiggle line is the response func- 
tion G^, wiggly line is the (p field and the vertex is the 
Si nonlinearity. The first diagram, is proportional to T, 
and is the graphical correction to Ai and T7, while the 
second one, survives even at zero temperature and renor- 
malizes 7. 

Expanding to quadratic order in the short-scale fields 
(j)'' and (f)'' and averaging we obtain 



^i>> 
Ai 
2 



[go(0<t-<t')]((</'r>t-0>t')'> 



r,t,t' 

+ 2go0<tSin[go(0<t " (brJ')Mr,t4>r,t')-^ 

« {l-qlC>{r = O,t = 0))SS<,<f'<] 



integrals over 6t (best evaluated in Fourier u spac e) and 
noting that the terms in last part of the Eq. g.58| renor- 
malize Ai, T7 and 7, respectively, we find 



dAY> = -Aiq^C> (r = 0,St = 0) 



-Ai 



Tq'o 



5(T7)W = -Aig2c>(r = 0,c. = 0) 



^7 



(1) - 



^2 -a n> 



Aig^za„G>(r = O,cj = 0) 
Ai^-^Ca-iK'-'dl 



(8.59) 



(8.60) 



(8.61) 



^2|^|3 

The calculation to second order in Ai given by 882 = 
— ^{SDy can of course be done directly, however, it is 
convenient to utili ze th e functional renormalization group 
calculation of Sec. V C . There we found for an arbitrary 
force correlation function A((/)) 



5So = 



■'r,t,t' 



x(A(<^<.-0rV)-A(O))^|i^d^ (S 

which when applied to the lowest harmonic A((p) 
Ai cos [got/)] gives 



.62) 



SS2 



r,t,t' 



cosh 



Cd-iA 



d-3 



4:K±"f\v 
(8.63) 



dl 



The diagram that leads to the above contribution to 6S2 , 
renormalizing the 0th and 1st harmonics, i.e. Aq, and 
Ai, respectively, is illustrated in Fig. 9 




^ Jv,t J St 



.t<~'tH'v,t 



C>(r = 0,<5t) 
5tG>{T = 0,5t) 



where in above we used causality (selecting the discretiza- 
tion with 9(0) = 0) and took advantage of the fact that 
the correlations functions are short-range in time to per- 
form small time gradient expansion. Performing above 



.57) ^^2 

Fig. 9: Diagram that contributes to the renormaliza- 
tion of A(u), to second order in A(u), i.e. renormalizes 
the 0th and 1st harmonics, Ag and Ai, respectively with 
the same notation as in Fig. 8. 

■ 58) therefore have to second order in Ai 

(2) _ A|g|Gd^-3 



5A, 



dl 



2„2/^ A<i-3 



(2) _ _ Afg^Oi-iA 



-dl 



(8.64) 
(8.65) 
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where Aq is the zeroth harmonic of the force-force cor- 
relator, i.e. the 0- ind epen dent correlator of the random 
drag. Combining Eqs.S.65 with the zeroth (tr ivial dimen - 



siona l rescaling) and first order results of Eqs. 8.47b , B.5£ - 
B.61 and rewriting the flow equations for the dimension- 
less couplings Ao, Ai, and T defined by Eqs.S.48-b.49 



we obtain the RG flow equations quoted in the main text, 



Sec. VI, Eqs. 6.4-6.6 



IX. APPENDIX F 

In this appendix we discuss the hydrodynamic equa- 
tions for the driven smectic of lines (illustrated in Fig. 4) 
that may be obtained in a three-dimensional supercon- 
ductor {dt — 2, di = 1). This smectic of lines has qual- 
itatively new properties as compared to the smectic of 
point particles described in Section IV. As discussed by 
Marchetti and NelsonEJ, in an isotropic flux-line liquid 
the conserved variables associated with hydrodynamic 
modes are the density field p(r) and the two components 
of a tangent field density T(r) = {Tx,Ty), describing the 
instantaneous bending of the lines away from the direc- 
tion of the external field (z direction). Since flux lines 
cannot start or stop inside the sample, the density and 
the two components of the tilt field are not independent 
dynamical variables, but are related by a continuity equa- 
tion in the time- like variable z, 

dJp+Vfr = 0. (9.1) 

This is simply the condition of no magnetic monopoles. 
The line smectic retains some degree of periodicity along 
the transverse direction y. This broken symmetry is de- 
scribed by the layer displacement (j)y = <f>. In addition 
both the conserved density and one component of the 
conserved tilt density are associated with independent 
hydrodynamic modes since they are not slaved to the 
layer displacement field. Having identified the relevant 
hydrodynamic variables for the line smectic as a one- 
dimensional layer displace men t (p, a density p and a tilt 
density t, related by Eq. 9.1, we now proceed to con- 



struct the phenomenological hydrodynamic equations for 
the line smectic and to study the spectrum of the hydro- 
dynamic modes of this system. 

The continuum hydrodynamic free energy for the over- 
damped line smectic is given by 



(9.2) 

where 5p = p — po, with po the equilibrium density. 
Here cl and C44 are the smectic bulk and tilt moduli, 
respectively, c\i is the in-layer compressibility, Kf and 
Ki are layer bending stiffnesses. The coupling constant 
K2 has dimensions of an elastic constant. The hydro- 
dynamic equations of the driven smectic contain addi- 
tional nonequilibrium terms, as compared to their equi- 
librium counterpart. The nonequilibrium terms can be 



constructed by preserving two important symmetries of 
the driven system, the invariance under inversions about 
the direction of the external drive {y —y, 4> —(f) 
and the broken translational invariance in the y direc- 
tions {y ^ y + a, (j) ^ (j) + a). 

Density and tilt density conservation requires the [den- 
sity and the tilt field to satisfy continuity equation£3. 



dtSp + V-j = 0, 



(9.3) 
(9.4) 



where j is the number current density and Jy is the an- 
tisymmetric tilt flux tensor. The density and tilt density 
fields are also related by the "continuity" equation 9.1. 



The equation for the layer displacement has the s ame 
structure as that for the two-dimensional lattice, Eq. L2 , 
and is repeated here for completeness 



{dt+vdx)(p = -2- — , 

P Po o(p 



(9.5) 



with To a kinetic coefficient. We will not discuss here the 
role of disorder on the hydrodynamics of the line smec- 
tic. Therefore we have not included any pinning force 
in the equations of motion. The hydrodynamic equa- 
tions need to be supplemented by constitutive relations 
for the current flux j and the tilt flux Jy . For simplicity 
we only consider here local hydrodynamics, but the non- 
locality of the elastic constants that is often important in 
flux-line systems can be trivially incorporated. The con- 
stitutive equations for the driven line smectic contain, 
however, new nonequilibrium terms not present in their 
equilibrium counterpart discussed in Ref. The two 
components of the current density are given by 



jx = (w + vi)6p + poV2dy(j) - aidzTx 



5Tis 
Sp 



jy = P0V3dx4> - a2dzTy - PqT2 dy 



5Tis 
Stx 

5p 



The antisymmetric tilt flux tensor is written as 



Ji'j — 



Vii ■ (vo X t) -f poVrT. • Vf X 



5t 



(9.6) 
(9.7) 

(9.8) 



All the parameters Vi and ai entering the nonequilibrium 
terms are proportional to the mean velocity v. Since the 
longitudinal part of the tilt vector can be eliminated in 
favor of the density using Eq. 9.1, it is convenient to work 
in Fourier space. We introduce longitudinal and trans- 
verse components of the tilt vector as 



T(q) = qtr/(q) -f (z x qt)rT(q), 



(9.9) 



with qt = qt/g*. Then n, = q* • t and tt = (z x q^) • t. 
By in sert in g th e constitutive equation for the fluxes in 
Eqs. dJ), (O) and (U), we obtain. 
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dt - iviqx + Diql + D2ql - (g^-De + qlDj)ql Sp = PQV2qxqy(j) + (-De - Dj)qxqyqzqtTT, 



dt - ivsqx + Dsql + D^ql + D^ql 



-iqy{D5 



2D7, 



Sp 
Po 



iq^qx — TT 
po 



(9.10) 
(9.11) 



dt + iViqx + Dgql - {Deql + Drq^qf tt = Poqzqt{v2ql - Vs,ql)(i) 

+i(ui +V4)qyqzSp + ql{D(i - D'j)qxqy—5p. 

qt 



(9.12) 



Finally, the longitudinal part of the tilt density is sim- 
ply related to the density, 



TL = op 

qt 

The "velocities" vi, V2, V3 have been defined as 



V2 

V3 



V + Vi, 
V2 + W3, 
V-V3 



(9.13) 



(9.14a) 
(9.14b) 
(9.14c) 



The coefficients Di have dimensions of diffusion constants 
and are given by 



first mode corresponds to the permeation mode of smec- 
tic liquid crystals and describes the transport of mass 
across the layers that can occur in these systems with- 
out destroying the layer periodicity. The second mode 
describes long- wavelength deformations of the layers and 
governs the decay of displacement fluctuations. Finally, 
the third mode governs the relaxation of tilt fluctuations, 
that, like density, can occur both in and out of the layers, 
while preserving the line smectic periodicity. 



Di 


= TiCl/po, 


(9.15a) 


D2 


= ^2Cl/P0, 


(9.15b) 


D3 


= ToKi/pa, 


(9.15c) 


Di 


= [TqcIi - T2K2)/po, 


(9.15d) 


D5 


= {roK2 - T2Cl)/P0, 


(9.15e) 


De 


= ai - TiCii/po, 


(9.15f) 


D7 


= 02 - r2C44/po, 


(9.15g) 


Ds 


= roK^/po, 


(9.15h) 


D9 


= TrCii/PQ- 


(9.15i) 



By solving the hydrodynamic equations in the long 
wavelength limit, we can find the hydrodynamic eigen- 
frequencies that govern the relaxation of density, tilt and 
displacement fluctuations. All the modes are propagat- 
ing at finite velocities and are given by 
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For stability, in addition to the conditions stated for 
the two-dimensional smectic, we must have Dq < 0, 
Dh - > 0, Dg > and Dq + < 0. The 
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